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We provide an approximate analytical expression of the mass-radius relation of a Newtonian self- 
gravitating Bose-Einstein condensate (BEC) with short-range interactions described by the Gross- 
Pitaevskii-Poisson system. These equations model astrophysical objects such as boson stars and, 
presumably, dark matter galactic halos. Our study connects the non-interacting case studied by 
Rufhni & Bonazzola (1969) to the Thomas- Fermi limit studied by Bohmer & Harko (2007). For re- 
pulsive short-range interactions (positive scattering lengths), there exists configurations of arbitrary 
mass but their radius is always larger than a minimum value. For attractive short-range interactions 
(negative scattering lengths), equilibrium configurations only exist below a maximum mass. Above 
that mass, the system is expected to collapse and form a black hole. We also study the radius versus 
scattering length relation for a given mass. We find that equilibrium configurations only exist above 
a (negative) minimum scattering length. Our approximate analytical solution, based on a Gaussian 
ansatz, provides a very good agreement with the exact solution obtained by numerically solving a 
nonlinear differential equation representing hydrostatic equilibrium. Our analytical treatment is, 
however, easier to handle and permits to study the stability problem, and derive an expression of 
the pulsation period, by developing an analogy with a simple mechanical problem. 



I. INTRODUCTION: A BRIEF REVIEW 



Quantum mechanics predicts that the particles can be 
classified in two groups: fermions [11,0] and bosons [1, 
In early works, these systems were studied in the absence 
of interaction and fundamental concepts such as the Pauli 
exclusion principle for fermions and the Bose-Einstein 
condensation for bosons were evidenced. Fermi statistics 
forces the particles to occupy successively higher energy 
levels while bosons at zero temperature all condense into 
the ground state. In these ideal systems, the distribu- 
tion of particles is spatially homogeneous. However, for 
very massive systems, as in astrophysical situations, one 
must take self-gravity into account. In that case, the 
system becomes spatially inhomogeneous and centrally 
condensed. This leads to the concepts of fermion [H, Q 
and boson (rl-fioj stars. 

Soon after the discovery of the quantum statistics for 
fermions by Fermi [l[ (1926) and Dirac (1926), Fowler 
pT| (1926), in a paper entitled "Dense Matter", under- 
stood that white dwarf stars owe their stability t o th e 
quantum pressure of the degenerate electron gas 134 1. 
Therefore, gravita tion al collapse is avoided by Pauli's ex- 
clusion principle jl35l |. Fowler modeled a white dwarf 
star at zero temperature by a completely degenerate 
Fermi gas in gravitational equilibrium. He noted that the 
star is "strictly analogo us to one gigantic molecule in its 
lowest quantum state" |l36l |. As shown by Milne (1930) 
and Chandrasekhar (1931a) [l3|, the resulting struc- 
ture is equivalent to a poly trope of index n = 3/2 so that 
the mass-radius relation of classical white dwarf stars is 
Mi?3 = 91.9h^/mlG^{fiH)^ where me is the mass of the 
electron, H the mass of the proton and fi the molecular 
weight. Chandrasekhar (1931b,1935) [23,111 took special 
relativity into account and showed that the star becomes 
equivalent to a polytrope of index n = 3 in the ultra- 



relativistic limit. From that result, he deduced the exis- 
tence of a maximum mass Mch = 3.09(?ic/G)'^/^(^.ff)~^ 
above which there is no hydrostatic equilibrium [l37{ . 
This mass Mch = 5.75Mq/ is of the o rder of the solar 
mass {Mch = I.UMq for fi ^ 2) [T38| . Massive stars 
cannot pass into the white dwarf stage. They undergo 
gravitational collapse and become neutron stars in which 
gravitational contractio n is arrested by the quantum 
pressure of the neutrons [l39l | . The structure of neutron 
stars was studied by Oppenheimer & Volkoff (1939) [s^l 
using general relativity. They obtained a maximum mass 
M = 0.376{hc/ G)^^^m-^ = 0.7Mq (and a corresponding 
radius R = 9.6km and density p = 5 lO^^gcm^"^) above 
which no equilibrium state exists (the Chandrasekhar cal- 
culations, neglecting general relativity, lead to a larger 
value of the maximum mass AI = 5.75Mq). In that 
case, nothing prevents the gravitational contraction of 
the star which becomes a black hole. The mass-radius 
relation M{R) of neutron stars makes a spiral while the 
mass-central density relation M{po) presents damped os- 
cillations. The series of equilibria becomes unstable after 
the first mass peak which corresponds to the maximum 
mass [33 - [35j . The existence of a limiting mass for self- 
gravitating fermions, the Chandrasekhar mass, is striking 
because it can be expressed in terms of fundamental con- 
stants similarly to Bohr's radius of the atom [3^. This 
limiting mass arises because of relativistic effects. For a 
non-relativistic fermion star, we have the mass-radius re- 
lation MR^ = 91.9?f/m^G^ that is vahd for small masses 
M <^ Mch- For general relativistic fermion stars [s^], the 
maximum mass and the minimum radius can be written 
Mch = 0.376Af|,/m2 and R = 3.52{Mp/mflp where 
Mp = (hcjGfl'^ is the Planck mass and Ip ^ {fiG/c^f/'^ 
is the Planck length. The radius can be expressed 
in terms of the Compton wavelength of the fermions 
Ac = h/mc as i? = 3.52(Mp/to)Ac. For m - ICeV/c^, 
we get Mch ^ lO^^kg ~ M© and R ^ km. 
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The preceding results are valid at T = 0. The self- 
gravitating Fermi gas at non-zero temperature has been 
studied in the context of statistical mechanics by Her- 
tel & Thirring (1971) [13, Messer (1981) ^ and more 
recently by Chavanis (2002) [s^ (an exhaustive descrip- 
tion of the phase diagrams in canonical and microcanon- 
ical ensembles is presented in the review [40|). In these 
studies, the system must be enclosed within an artifi- 
cial "box" so as to avoid its complete evaporation. For 
large systems, one recovers the classical isothermal self- 
gravitating gas that may undergo a gravothermal 
catastrophe j4l| in the microcanonical ensemble (fixed 
energy) or an isothermal collapse in the canoni- 
cal ensemble (fixed temperature). If the particles are 
fermions, the gravitational collapse stops when quantum 
degeneracy effects (Pauli exclusion principle) come into 
play. This leads to the formation of a composite struc- 
ture made of a completely degenerate and very compact 
nucleus (fermion ball) surrounded by a dilute atmosphere 
jiof . Self-gravitating fermions were also discussed in cos- 
mology, in models where dark matter is made of mas- 
sive neutrinos [l40] . Originally, the self-gravitating Fermi 
gas with neutrino masses in the ~ cV/c^ range was pro- 
posed by Ruffini and co-workers (1982-1990) [44| - |4a | as a 
model for dark matter halos (i? ~ 10^'^cm ~ lOOkpc and 
M ^ W^^g ^ IO^^Mq) and clusters of galaxies. Then, 
Viollier and coworkers (1997-2001) [i^-Iii suggested that 
degenerate superstars composed of weakly interacting 
fermions in the ^ lOkeV/c^ range could be an alterna- 
tive to the supermassive black holes that arc reported to 
exist at the centers of galaxies (e.g. M = 2.6 10^ Mq and 
R = 18mpc in our Galaxy). Finally, Bilic et al. (2003) 
jsof showed that a weakly interacting fermionic gas at fi- 
nite temperature could provide a self-consistent model of 
dark matter that describes both the center and the halo 
of the galaxies. Since the density of a self-gravitating 
isothermal sphere decreases as at large distances, this 
model is consistent with the flat rotation curves of the 
galax;ies. On the other hand, since the core is degenerate 
in the sense of quantum mechanics (Pauli exclusion prin- 
ciple) , it leads to flat density profiles and avoids the cusp 
problem of cold dark matter (CDM) models (see below). 
In addition, the gravitational collapse of fermionic mat- 
ter leads to a compact object (fermion ball) at the center 
of the galaxy that could mimic a central black hole. 

Coincidentally, the Fermi-Dirac statistics also arises in 
the theory of violent relaxation developed by Lynden-Bell 
(1967) l5l| for the Vlasov-Poisson system. However, in 
that case, the origin of the "degeneracy" is due to dynam- 
ical constraints (Liouville's theorem) instead of quantum 
mechanics (Pauli's principle). This theory was initially 
developed to describe collisionless stellar systems such as 
elliptical galaxies. In that case, the non-degenerate limit 
may be the most relevant [5l|. However, as proposed 
by Kull et al. (1996) [S^ and Chavanis & Sommeria 
(1998) [l^l, this approach (with dynamical degeneracy 
retained) could also apply to dark matter halos and pro- 
vide a much more efficient relaxation mechanism than in 



the fermion scenario. Indeed, the violent relaxation of 
collisionless systems (leading to the Lynden-Bell statis- 
tics) takes place on a few dynamical times while the coUi- 
sional relaxation of fermions (leading to the Fermi-Dirac 
statistics) is larger than the age of the universe by sev- 
eral orders of magnitude. Therefore, it is not clear how 
the fermions have thermalized and how they can possess 
sufficiently large temperatures. By contrast, the Lynden- 
Bell theory predicts a large effective temperature (even 
if T = initially), a density profile at large distances 
consistent with the fiat rotation curves of galaxies and 
an effective exclusion principle at short distances that 
could avoid the cusp problem and lead to fermion balls 
mimicking black holes, just like in the fermionic scenario. 

Therefore, fermion stars rapidly found applications in 
relation to white dwarf stars, neutron stars, massive neu- 
trinos in dark matter models, and in the violent relax- 
ation of collisionless self-gravitating systems. 

The concept of boson stars was introduced by Kaup 
(1968) [sl and Ruffini & Bonazzola (1969) ^ ahhough 
no astrophysical application of these objects was known 
at that time. They were just hypothetical stars whose 
main interest was ruled by the fundamental laws of 
physics that govern their structure. In a sense, boson 
stars are descendent of the so-called geons of Wheeler 
(1955) [13 except that they are built from scalar parti- 
cles (spin-0) instead of electromagnetic fields, i.e. spin-1 
bosons. Kaup and Ruffini & Bonazzola considered the 
T = limit jl4ll | at which the bosons form a Bose- 
Einstcin condensate (BEG) and showed that the con- 
cept of an equation of state breaks down. Indeed, at 
zero temperature a boson gas has a vanishing pressure 
and is unable to sustain any configuration of equilibrium 
in the classical perfect fiuid approximation. But if one 
considers a quantum system of massive self-gravitating 
bosons in their ground state self-consistently, stable equi- 
librium configurations exist. In that case, all the bosons 
are in the same quantum state and they are described by 
a unique wave function ■0(r). Therefore, boson stars can 
be regarded as macroscopic quantum states that arc only 
prevented from collapsing gravitationally by the Hcisen- 
berg uncertainty principle AxAp ~ h. In the Newto- 
nian approximation, a self-gravitating BEG is described 
by the Schrodinger-Poisson system and in the relativistic 
case one must couple the Klein- Gordon equation to the 
Einstein field equations. The Newtonian approximation 
is valid for sufficiently small masses and yields the mass- 
radius relation MR = 9.9h^ /Gm? (where R is the radius 
containing 99% of the mass) [s^. The radius decreases 
as mass increases, like for classical white dwarf stars, 
but the scaling is different. When relativistic effects are 
taken into account, there exists a maximum mass, the 
Kaup mass Mxaup = 0.633Afp/?Tt, above which no equi- 
librium configuration exists [54l |55|. In that case, the 
system collapses to a black hole. The minimum radius 
Rmin — 6.03Ac corresponding to the Kaup mass is of 
the order of the Gompton wavelength Ac = fi/mc of the 
bosons. These results were re-derived independently by 
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Thirring (1983) [H] and Breit et al. (1984) 

There exists remarkable similarities between boson and 
neutron stars [U, iKl, .6i|. For example, the mass- 
central density relation M{pq) of boson stars exhibits 
damped oscillations and the series of equilibria becomes 
dynamically unstable after the first mass peak corre- 
sponding to the maximum mass [g^ - Ib^ ]. On the other 
hand, the mass-radius relation M{R) has a snail-like 
(spiral) structure, and the mass-particle number rela- 
tion M{N) presents cusps. The Newtonian approxima- 
tion is valid for sufficiently low densities corresponding 
to M <C MKaup and R ^ Xc- On the other hand, 
boson stars become relativistic when they approach the 
maximum mass Mxaup = 0.633Afp/m and the corre- 
sponding radius Rmin = 6.03Ac. Except for matters of 
scales, all these results are remarkably similar to those 
obtained for neutron stars [33 - l35j . However, there also 
exists crucial differences between self-gravitating objects 
made of fermions or bosons. In particular, boson stars 
are stopped from collapsing by Heiscnberg's uncertainty 
principle while for fermion stars gravitational collapse is 
avoided by Pauli's exclusion principle. This difference 
is reflected in the critical mass of stable configurations: 
Mch Mp/m? for fermions and Mxaup Mp/m for 
bosons. This difference is due principally to the fact 
that all the N bosons are in the ground state whereas 
the fermions, according to the Pauli exclusion prin- 
ciple, are distributed in the N lowest energy states of 
the phase space. As a result, the mass of boson stars 
MKaup ~ Mp /m is generally much smaller than the mass 
of fermion stars Mch ~ Mp/iri^. They differ by a factor 
m/Mp <^ 1. For example, for m ~ IGeV/c^ for which 
m/Mp ^ 10~^^, one can estimate the total mass of a 
boson star to be M ^ lO^^kg ^ 10~^^AfQ and its radius 
R ~ 10~^^m yielding a density 10'^* times that of a neu- 
tron star. Contrary to the mass of neutron stars that is of 
the order of a solar mass, the mass of these mini boson 
stars is too small to be astrophysically relevant. They 
could play a role, however, if they exist in the universe in 
large quantity or if the mass m of the bosons is extraor- 
dinary small leading to macroscopic objects with a mass 
comparable to the mass of the sun (or even larger) [lo| . 
This is the case, in particular, for axionic boson stars that 
could account for the mass of MACHOs (between 0.3 and 
0.8 Mq) if the axions have a mass m ~ 10^^°eV/c^ [65j . 

Colpi et al. (1986) [g^ considered the case of self- 
interacting scalar fields with a iA|(/)|^ term and found 
that the resulting configurations differ markedly from the 
non-interacting case even when A 1. In that case, the 
maximum mass may be comparable 

with the Chandrasekhar mass of self-gravitating fermions 
when A ~ 1. Similarly, the radius R ^ ^f\{Mp /m)\c of 
self-interacting boson stars may be much larger than the 
Compton wavelength and become comparable with the 
radius of fermion stars. For m ^ IGeV/c^, boson star 
configurations exist with a mass M ~ lO^^VAkg and a 
radius R ^ -y/Akm similar to those of neutron stars. For 
smaller masses m ^ IMcV/c^, we get M - 1036\/Akg 



and R ^ lO^V^km. The radius is comparable to that 
of the sun but it encloses 10^ solar masses. These pa- 
rameters are reminiscent of supermassive black holes in 
Active Galactic Nuclei, so that boson stars, like fermion 
stars, could be an alternative to black holes [G^I- These 
enhancements of mass and size are due to the parameter 
\f\Mplm that is very large even for small A due to the 
smallness of m relative to Mp. Therefore, self-coupling 
can significantly change the physical dimensions of bo- 
son stars, making them much more astrophysically inter- 
esting. The self-interaction has the same effect on the 
bosons as the exclusion principle on fermions. It plays 
the role of an interparticle repulsion (for A > 0) that 
dominates over uncertainty pressure and prevents catas- 
trophic gravitational collapse. Colpi et al. [HBl showed 
that a quartic self-interaction is equivalent to an effective 
barotropic pressure. At low densities p = (A?i"^/4m'*c^)e^ 
(polytropc n = 1) and at high densities p = e/3 (linear) 
like in the core of neutron stars [33 - !35l | (this result is not 
explicitly given in (66j but it can be obtained as a limit 
of their Eq. (17)). This strengthens the analogy between 
boson stars and neutron stars. Since boson stars may 
have masses comparable to the mass of neutron stars, or 
even larger, they become astrophysical relevant and may 
play a role in the problem of dark matter [g^ [g^ . 

A wide "zoology" of exotic particles that could form 
dark matter has been proposed. In particular, many 
grand unified theories in particle physics predict the ex- 
istence of various exotic bosons (e.g. axions, scalar neu- 
trinos, neutralinos) that should be present in consider- 
able abundance in the universe and comprise (part of) 
the cosmological missing mass [H, [g^. Although the 
bosonic particles have never been detected in accelerator 
experiments, they are considered as leading candidates 
of dark matter and might play a significant role in the 
evolution and the structure of the universe. The forma- 
tion of boson stars was investigated by Madsen & Liddle 
(1990) [t^ and is now relatively well understood. A spa- 
tially homogeneous distribution of self-gravitating bosons 
can undergo a sort of Jeans instability as describe d by 
Khlopov et al. (1985) [7l| and Bianchi et al. (1990) [zf. 
When the perturbation has sufficiently grown, the cloud 
collapses under its own gravity at first in free fall. Then, 
as nonlinear gravitational effects become important at 
higher densities, the configuration starts to oscillate and 
settles into a compact bosonic object through the radia- 
tion of the scalar field. This gravitational cooling process 
has been evidenced and studied by Seidel & Suen (1994) 
[t^. This is a dissipationless mechanism similar in some 
respect to the violent relaxation of collisionless stellar 
systems [5l|, [tI] but ending on a unique final state (bo- 
son star) independent on the initial conditions. There- 
fore, (mini) boson stars could be the constituents of dark 
matter halos [g^, [g^. Later, it was suggested that dark 
matter halos themselves could be gigantic self-gravitating 
BECs. This idea was introduced in order to solve basic 
problems inherent to cold dark matter models. 

Dark matter is one of the most important puzzles 
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in modern physics and cosmology. Cold dark matter 
(ACDM) models with a small cosmological constant com- 
prising weakly interacting massive particles (WIMPs), 
such as the lightest neutralinos, are presently favored by 
theorists over hot dark matter (HDM) models compris- 
ing relativistic light neutrinos. Although the cold dark 
matter (CDM) model is popular and remarkably success- 
ful in explaining the large-scale structure of the universe 
(75I , it seems to encounter many problems on the scale of 
galactic or sub-galactic structures. Indeed, CDM simu- 
lations lead to cuspy density profiles at galactic cen- 
ters (in the scales of the order of 1 kpc and smaller) [t^ 
while most rotation curves indicate a smooth core density 
[77I . On the other hand, the predicted number of satellite 
galaxies around each galactic halo is far beyond what we 
sec around the Milky Way [tI]. These problems might 
be solved, without altering the virtues of CDM models, if 
the dark matter is composed of scalar particles in a cold 
BEC. The wave properties of the dark matter may stabi- 
lize the system against gravitational collapse as a conse- 
quence of the Heisenberg uncertainty principle, providing 
halo cores instead of cuspy profiles. The resulting coher- 
ent configuration may be understood as the ground state 
of some gigantic bosonic atom where the ultra- light boson 
particles are condensed in a single macroscopic quantum 
state V'(r). In these models, the formation of dark matter 
structures at small scales is suppressed by the quantum 
uncertainty principle. This property could alleviate the 
problems of the CDM model such as the cusp problem 
and the missing satellite problem. 

The first suggestion that galactic halos are formed by 
bosons, either in their quantum ground state (BEC) or in 
an appropriate isothermal distribution, is due to Baldc- 
schi et al. (1983) [t^. In order to yield masses and sizes 
that agree with those of galactic halos (M ^ 10*^ g and 
R ~ 10^3 cm), the mass of the bosons must be extremely 
small m ~ 10^^"'eV/c^ (estimated with the Newtonian 
mass-radius relation) . The idea that galactic halos could 
be a giant system of "Bose liquid" was also advanced 
by Sin (1994) (sOj who studied rotation curves induced 
by self-gravitating BECs. Sin considered the ultralight 
pseudo Nambu-Goldstonc boson (axion) appearing in the 
late-time cosmological phase transition theories as a ma- 
jor dark matter candidate. Since it is almost massless 
m ~ 10~^^eV/c^, its nature is more wavelike than par- 
ticlelike so that the dark matter distribution must be 
treated quantum mechanically. Therefore, galactic halos 
can be considered as self-gravitating Bose liquids whose 
collapse is prevented by the uncertainty principle. At the 
galactic scale, the Newtonian approximation turns out to 
be relevant. Sin considered condensation wave functions 
for the galactic halo that have nodes (excited states) be- 
cause the rotation curve of the zero node solution falls 
too fast to explain the flatness of the rotation curve of 
many galaxies. Because of the systems' wavelike nature, 
he found that the rotation curve of galaxies has a ripplc- 
likc fine structure that seems to agree with observations. 
However, this may be a coincidence because the presence 



of wiggles in the rotation curves is rather related to spi- 
ral arms inside the disks, and excited states are generally 
unstable ■ Indeed, the excited modes decay to the 
ground state through emission of gravitational radiation, 
a process similar to atomic transitions [8l|, [s^ . In this 
sense, boson stars arc like gravitational atoms. 

An alternative scalar field matter model for dark ha- 
los of galaxies was developed by Schunck (1998) [83| . 
He obtains rotation curves that give a very good agree- 
ment with observations of spiral and dwarf galaxies. In 
this model, the density decreases at large distances like 
yielding asymptotically flat rotation curves. Fur- 
thermore, this model produces oscillations around this 
asymptotic value that can match, in some cases, the data. 
Since these oscillations ( "wiggles" ) do not correspond to 
excited states, the solutions are stable. Although these 
results are obtained in the Newtonian limit, it is shown 
that the radial pressure coming from general relativity 
plays a role for the rotation velocity that is comparable 
to the contribution of the normal part. On the other 
hand, in the strongly relativistic case, this model leads 
to massive objects with large redshift values and rotation 
velocities. This might be an explanation for the large en- 
ergy contributions seen in quasars |83|. At about the 
same period, Matos & Guzman (1999) [SJl introduced a 
model in which dark matter is a scalar field embedded 
in an exponential (and later cosh) scalar potential. This 
model produces a density profile of the form l/(r^ + ^c) 
which gives a good agreement with the rotation curve of 
galaxies and accounts for their flatness at large distances 
j8a | . It also explains the suppression of subgalactic struc- 
tures (since it produces a sharp cut-off in the mass power 
spectrum) and the smoothness of galaxy core halos [s^ . 
This model has only one free parameter, the scalar field 
mass, whose determined value m ~ 10~^'^eV/c^ accounts 
both for the typical mass of galactic halos and for cos- 
mological observations (two a priori independent mea- 
surments) ISTll . A similar mass was found by Arbey et 
al. (2001) jig and Silverman & Mallett (2002) 

In these models, the self-interaction of the particles is 
neglected and the mass of the bosons must be extremely 
small in order to reproduce the characteristic mass and 
size of galactic halos. Such an ultralight scalar field (e.g. 
an axion) with m ~ 10~^^eV/c^ was called "fuzzy cold 
dark matter" (FCDM) by Hu et al. (2000) ^ who dis- 
cussed its overall cosmological behavior. An alternative 
to this unnatural small mass is to take self-coupling into 
account. Indeed, in the context of boson stars, the work 
of Colpi et al. (1986) |6^ has demonstrated that, for 
the same value of the boson mass m, even a small cou- 
pling can considerably change the mass and size of self- 
gravitating BECs. In that case, values of the mass in 
the ballpark of an eV/c^ may be compatible with a size 
of a few kiloparsecs. Lee & Koh (1996) [9lj investigated 
relativistic boson stars with a self-interacting scalar field 
as a model of galactic halos. A massive scalar field or 
a boson condensate with quartic - or close to quartic - 
self-coupling was also proposed as a possible dark mat- 
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ter candidate by Peebles (2000) who called it "fluid 
dark matter" and by Goodman (2000) [1^ who called it 
"repulsive dark matter". Similarly, Arbey et al. (2003) 
[q^ I considered a self-coupled charged scalar field which 
is equivalent to a self-gravitating Bose condensate. For 
m^/X ^ 50 (eV/c^)''^, they obtained a very good agree- 
ment with the measurements of the circular speed of the 
dwarf spiral DD0154. Bohmer & Harko (2007) [H pur- 
sued the idea that dark matter is in the form of a self- 
gravitating BEG and studied the condensate by using 
the non-rclativistic Gross-Pitaevskii equation coupled to 
the Poisson equation. They took self-interaction into ac- 
count via a quartic nonlinearity and used the Thomas- 
Fermi approximation which becomes exact for N — > +oo. 
Under these assumptions, the BEG is equivalent to a 
barotropic fluid with a polytropic equation of state of 
index n = 1 [9ll, [93l - l95j . This leads to a length-scale 
R = Triafi'^ /Gm?)^/"^ for bound objects that is indepen- 
dent on their mass (a is the scattering length). Bohmer 
& Harko considered the example of a galactic dark mat- 
ter halo extending up to i? = 10 kpc = 3.08 10^^ cm with 
a mass of the order of M = 3 10^^ Mq yielding an aver- 
age density p = 5.30 10"^"' g cm~^. For a = 5.7710~^cm 
(a typical value in terrestrial BEG experiments (9^) 
they found that the mass of the particles forming the 
condensate dark matter halo is of the order of the eV 
(more precisely, m = 1.44 eV/c^ for a = 10^ fm and 
m = 14meV/c^ for a = 1 fm). They also determined the 
rotation curves created by a BEG and found a very good 
agreement with the observational data for several low 
surface brightness galaxies. For recent studies of scalar 
field/BEG dark matter models see, e.g., Refs. [97l - [l04 |. 

This detailed introduction shows that the idea of bo- 
son stars has a long and rich history and that they may 
play a role in different areas of astrophysics. If we con- 
sider dark matter halos, the Newtonian approximation is 
sufficient. In previous works, two limits have been consid- 
ered. Ruffini & Bonazzola [5al neglect the self-interaction 
of the particles (a = 0) and solve the Schrodinger-Poisson 
system. In that case, the equilibrium state results from 
the balance between the gravitational attraction and the 
Heisenberg principle equivalent to a quantum pressure. 
Alternatively, Bohmer & Harko (9^ take into account 
the self-interaction of the particles (a > 0) and study 
the Gross-Pitaevskii-Poisson (GPP) system. They con- 
sider the Thomas-Fermi (TF) limit which amounts to 
neglecting the quantum pressure. In that case, the equi- 
librium state results from the balance between the grav- 
itational attraction and the small-scale repulsion due to 
scattering. As we shall see, the TF limit is valid if 
GN^m^a/h^ ^1. In the present work, we shall connect 
these two limits by considering the general case where 
both short-range interactions and quantum pressure are 
taken into account. We shall also treat the case where the 
self-interaction is attractive instead of repulsive. Since 
atoms may have negative scattering lengths in terrestrial 
BEG experiments [96|, it may be useful to consider the 
possibility of attractive interactions in our general study. 



The TF approximation cannot be employed in that case 
and we have to use the complete set of equations. At- 
tractive self-interaction is equivalent to a negative pres- 
sure p = — that adds to the gravitational attraction. 
In that case, we find the existence of a maximum mass 
Mmax = 1.012h/ y^\a\Gm above which the system can- 
not be in equilibrium. This maximum mass can be very 
small, as small as the Planck mass Mp (!), meaning that 
when the self-interaction is attractive the system is very 
unstable. When applied to a cosmological context jl06l |. 
an attractive self-interaction could enhance the gravita- 
tional collapse and accelerate the formation of structures 
in the universe. 

The paper is organized as follows. In Sees. HI] and 
mil we provide general results concerning the Gross- 
Pitaevskii-Poisson system. We specifically consider the 
non-interacting case and the Thomas-Fermi limit. In 
Sec. IIVI we obtain an analytical approximate expression 
of the mass-radius relation of self-gravitating BEGs with 
positive or negative scattering lengths by using a Gaus- 
sian ansatz for the wave function and developing a simple 
mechanical analogy. In Sec. |Vl wc study the Jeans in- 
stability of an infinite homogeneous self-gravitating BEG 
by taking into account the self- interaction of the part icles 
that was ignored in previous works. In paper II jl05l |. we 
show that our approximate analytical approach gives a 
good agreement with the exact results obtained by nu- 
merically solving the equation of hydrostatic equilibrium. 
In Paper HI, we extend our analytical approach to more 
general situations. 

II. THE GROSS-PITAEVSKII-POISSON 
SYSTEM 

A. The mean-field Gross-Pitaevskii equation 

Following Bohmer & Harko [1^ , we model dark matter 
halos as a self-gravitating Bose-Einstein condensate with 
short-range interactions. Since the cosmic BEG has a rel- 
atively low mean mass density, we can use the Newtonian 
appr oximation. At T = 0, all the bosons have condensed 
[l42l | and the system is described by one order parameter 
ip{r, t) called the condensate wave function. In the mean- 
field approximation, the ground state properties of the 
cond ensa t e ar e described by the Gross-Pitaevskii equa- 
tion flOTlflOSt : 

dib fi^ 
ih^{v,t) = - — Ai;{r,t) +m^tot{r,t),p{r,t), (1) 

*tot(r,t) = J p{r',t)u{\r~r'\)dr\ (2) 
p(r,t) = 7Vm|^(r,t)|2, (3) 
I mr,t)\' dr ^ 1. (4) 
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Equation (j4]) is the normalization condition, Eq. ([3|) gives 
the density of the BEC, Eq. ([2]) determines the associ- 
ated potential and Eq. ([T|) determines the wave function. 
We assume that the potential of interaction can be writ- 
ten as u = ulu + usb. where u^n refers to the long-range 
gravitational interaction and usr to the short-range in- 
teraction. We assume that the short-range interaction 
corresponds to binary collisions that can be modeled by 
the pair contact potential usR{r—r') ~ gS{r~r'), where 
the coupling constant (or pseudo-potential) g is related 
to the s-wave scattering length a through g = ATrah^ /m^ 
[96j . For the sake of generality, we allow a to be positive 
or negative (a > corresponds to short-range repulsion 
and a < corresponds to short-range attraction). Un- 
der these conditions, the total potential can be written 
^tot = ^ + h{p) where $ is the gravitational potential 
and 



Kp) = gp = gNm\ip\'^, 



(5) 



is an effective potential modeling short-range interac- 
tions. When this form of potential is substituted in Eq. 
(21), we obtain the Gross-Pitaevskii-Poisson (GPP) sys- 
tem 



in-— = Ayj 

ot 2m 



A$ = AnGp = 47rG7Vm|?/>p 



(6) 



(7) 



that will be the object of focus in this paper. In the gen- 
eral formalism developed in the sequel, we will consider 
an arbitrary potential h{p). However, for specific appli- 
cations, we will consider the potential ([S]). More general 
situations will be studied in Paper III. 



B. The Madelung transformation 

Let us use the Madelung |l09l | transformation to 
rewrite the GPP system in the form of hydrodynamic 
equations. From the wave function 



(8) 



where A{r, t) and S{r, t) are real functions, we introduce 
the density and velocity fields 

p^ NmlM^ ^ NmA^, u=—VS. (9) 

m 

We note that the flow defined in this way is irrotational 
since V x u = 0. Substituting Eq. ([8]) in Eq. ([6]) and 
separating real and imaginary parts, we obtain 



dS 



1 



dt 2m 



|+V.(,u) = 0, 



{VSf + m<^ + mh{p) 



0, 



(10) 



(11) 



where 



Q 



2m y/p 



Am 



Ap 
P 



l{Vpl 

2 fl2 



(12) 



is the quantum potential. The first equation is similar to 
the equation of continuity in hydrodynamics. The second 
equation has a form similar to the classical Hamilton- 
Jacobi equation with an additional quantum term. It can 
also be interpreted as a generalized Bernouilli equation 
for a potential flow. Taking the gradient of Ec^. (fTTj) 
and using the well-known identity (u • V)u = V(u2/2) — 
u X (V X u) which reduces to (u • V)u = V{u'^/2) for 
an irrotational flow, we obtain an equation similar to the 
Euler equation with an additional quantum potential 



- + (u.V)u = 



-V/i- V4>- — VQ. (13) 
m 



This equation shows that the effective potential h appear- 
ing in the GP equation can be interpreted as an enthalpy 
in the hydrodynamic equations. We can rewrite Eq. (jl3p 
in the form 

Qw 1 1 

^ + (u-V)u = — Vp-V$ VQ, (14) 

at p m 

where p{r,t) is a pressure. Since h{r,t) = h[p{r,t)], the 
pressure p{r,t) = p[p{r,t)] is a function of the density 
(the flow is barotropic). The equation of state p{p) is 
determined by the potential h(p) through the relation 



h'{p) =. 



P'ip) 



(15) 



This yield s p{p) = ph{p) — H{p) where _ff is a primitive 
of h |143| . In conclusion, the GPP system is equivalent 
to the "hydrodynamic" equations (HUl), (Ull) and 0. We 
shall refer to these equations as the quantum barotropic 
Euler equations. In the classical limit ?i — 0, the quan- 
tum potential disappears a nd w e recover the ordinary 
barotropic Euler equations |llCll | . For a potential of the 
form ([S]), the equation of state is 



This is equivalent to a polytropic equation of state 



(16) 



P = Kp'^, 



1 

n 



(17) 



with a polytropic constant K = ^Tioh^ jw? and a poly- 
tropic index n = \ (i.e. 7 = 2). Inversely, the effective 
potential associated with the general polytropic equation 
of state (Ull) is h{p) = {K-f j ia - \)\p'^-^ . 

Remark: the quantum pote ntial (|12[) first appeared 
in the work of Madelung |l09l | and was rediscovered by 
Bohm jlll| (it is sometimes called "the Bohm poten- 
tial"). We note the identity 



-VQ 



-d P 



(18) 
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where Pij is the quantum stress (or pressure) tensor 



Am? 



pdidj Inp, 



or 



(19) 



(20) 



This shows that the quantum potential is equivalent to 
an anisotropic pressure. 



C. The time-independent GP equation 

If we consider a wave function of the form 

V(r,<) = A(r)e-'^^, (21) 
we obtain the time-independent GP equation 



the quantum potential and the repulsion (for a > 0) or 
the attraction (for a < 0) due to the short-range inter- 
action (scattering). This equation is equivalent to Eq. 
(1221) • Indeed, integrating Eq. (gSI) using Eq. ([T5|), we 
obtain Eq. pS]) where the eigenenergy E appears as a 
constant of integration. Combining Eq. (|25l) with the 
Poisson equation ([7]), we obtain the fundamental equa- 
tion of hydrostatic equilibrium with quantum effects 



4a (^] = A.Gp. 
2m^ \ ^ J 



(26) 



For an equation of state of the form (fT6|). it becomes 



—^P' 



2m 



V VP J 



(27) 



There are two important limits that we discuss in the 
following. 



- ^^A^{r) + m(<I>(r) + h{p))^{r) = E^{v), (22) 

where ip{r) = A{r) is real and p(r) = Nm'i}P'{Y). Dividing 
Eq. ^ by V'(r), we get 



m$ -|- mh{p) 



2m y/p 



or, equivalently, 



m<I> -I- mh{p) + Q = E. 



(23) 



(24) 



This relation can also be obtained from the quantum 
Hamilton- Jacobi equation ([TT|) by setting S = —Et. 
Combined with the Poisson equation ([7|), we obtain an 
eigenvalue equation for the wave function ?/'(r) where the 
eigenvalue E is the energy. In the following, we shall be 
interested by the fundamental eigenmode corresponding 
to the smallest value of E. For this mode, the wave func- 
tion ipir) is spherically symmetric and has no node so 
that the density profile decreases monotonically with the 
distance. 



D. Hydrostatic equilibrium 



The time- independent solution ([25)1 can also be ob- 
tained from the quantum barotropic Euler equation (jl4p 
since it is equivalent to the GP equation. The steady 
state of the quantum barotropic Euler equation (HH) , ob- 
tained by taking dt = and u = 0, satisfies 



(25) 



This generalizes the usual condition of hydrostatic equi- 
librium by incorporating the contribution of the quan- 
tum potential. Equation (|25p describes the balance be- 
tween the gravitational attraction, the repulsion due to 



E. The non-interacting case 

The non-interacting case corresponds to a = 0. This 
is the situation first considered by Ruf fini fc Bonazzola 
[ssf and revisited by Membrado et al. 112 with a dif- 



ferent method. In that case, the condition of hydrostatic 
equilibrium (P7)) reduces to 



2m? \ yfp J 



= ATlGp. 



(28) 



This corresponds to the balance between the gravita- 
tional attraction and the repulsion due to the quantum 
pressure arising from the Heisenberg uncertainty princi- 
ple. This equation can be solved numerically to yield the 
density profile [s^ Ill2| . The den sity decays smoothly to 
infinity (see, e.g.. Fig. 1 in jll2{ ) so that its support is 
not compact, contrary to fermion stars at T = (5[. The 
radius of the configur ation containing 99% of the mass 
has been computed in [ll2l |. They obtained the value 



i?99 = 9.9- 



GMm^ 



(29) 



This radius is much smaller than the gravitational Bohr 
radius as = t?/{Gm^) by a factor 1/7V ^ 1. They also 
found that the quantum kinetic energy Gg, the potential 
energy W, the total energy Efot and the eigenenergy E 
(see their definitions in Sec. IIII Ap are given by 



0.05426 



G^Afim^ 



W = -0.10852 



G^lVPm^ 



Et. 



-0.05426 



G'^M^m^ 



(30) 
(31) 
(32) 
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E = -0.16278- 



(33) 



Finally, they determined numerically the rotation curve 



Voir) 



GM{r) 



1/2 



(34) 



prod uced by a Newtonian self-gravitating BEC (see Fig. 
2 in |ll2| ) and mentioned application to galactic dark 
matt er halos. To our knowledge, Membrado et al. (1989) 
[ll2| were among the first authors in the literature to 
propose that dark matter halos could be a condensate 
boson sphere and to compute the corresponding rota- 
tion curve. However, they did not take into account self- 
coupling which is necessary to obtain physically relevant 
results. 



F. The Thomas-Fermi approximation 

The Thomas-Fermi (TF) approximation amounts to 
neglecting the quantum potential in Eq. ([25l) . This is the 
situation considered by Bohmer & Harko j9^ in relation 
to dark matter halos with repulsive self-interaction a > 0. 
In that case, Eq. reduces to the usual condition of 
hydrostatic equilibrium 



Vp + pV$ = 0. 



(35) 



This corresponds to the balance between the gravita- 
tional attraction and the repulsion due to the short-range 
interaction. Combined with the Poisson equation ([7]), we 
obtain the fundamental equation of hydrostatic equilib- 
rium 



V • ( ^ ) = 47rGp. 



(36) 



For an equation of state of the form (|16p , it can be rewrit- 
ten 



-p = 0. 



(37) 



This equation, which is equivalent to the Lane-Emden 
equation for a polytrope of index n = 1, can be solved 
analytically The density profile is given by the for- 
mula 



p(r) = sm ( — 

7rr \ n / 

where po is the central density and 

-2 \ 1/2 



i? = TT 



Gm? 



(38) 



(39) 



is the radius of the configuration at which the density 
vanishes (compact support). The radius of a polytrope 



n = 1 is independent on the mass M These re- 
sults have been derived by several authors in the context 
of self-gravitating BECs [al IH-iil using different for- 
malisms. The radius containing 99% of the mass is given 
by i?99 = 0.9542421 li?. The central density is deter- 
mined by the mass according to 



ttM 

IR3 



M 
4^ 



G^3x3/2 



ah-" 



(40) 



Using the analytical expression ([38)) of the density profile, 
we find that the moment of inertia (j72[) and the internal 
energy (|63|) are given by 



/ = 



1-A 

^2 



MR^ = (tt^ - 6) 



Gr 



U = 



GAP 



1 g^V/^aP 



An 



(41) 



(42) 



On the other hand, in the TF approximation, the steady 
state equation reduces to 



771$ 



-p = E. 



(43) 



It can be used to determine the eigenenergy E. Indeed, if 
we evaluate this relation at r = i? at which $ = ~GM /R 
and (0 = 0, we find that 



E = -- 



GMt 



R 



1 G^Prn^PM 



(44) 



Then, Eq. (|43]) with Eqs. (|38l) and dH]) determine the 
gravitational potential $(r). On the other hand, multi- 
plying Eq. (|43p by p, integrating over the configuration, 
and using the expressions (j42|) and (|44|) of the internal 
energy U and eigenenergy i?, we find that the potential 
energy (|64p is given by 



W 



3GA/2 



4i? 



3 G3/2„^3/2j^^2 

'4^ ai/2?i 



The total energy is 

Etot = U + W ^ 



GAP 
2R 



1 G^Pm^PAd^ 



2tt 



i^Pn 



(45) 



(46) 



Finally, using Eq. ((38)) . the rotation curve (|34p has the 
analytical expression 



AGpoR^ 



'R . 


/'TTr\ 


/7rr\ 


— sm 






nr 




-'°'{r)_ 



(47) 



For r — > 0, the velocity increases linearly with r as for a 
uniform sphere with density po: Vc{r) ^ {4:T:poG/3)^Pr. 
For r > R, we recover the Keplerian law Vc{r) = 
(GAf/r)^/^. The rotation curves created by dark matter 
halos made of self-gravitating BECs in the TF limit have 
been studied by Arbey et al. [g^l and Bohmer & Harko 
[95| who showed that they provide a good agreement with 
measured rotation curves of certain spiral galaxies. 
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G. Dimensional analysis and validity of the TF 
approximation 

In the absence of short-range interaction, the struc- 
ture of the self-gravitating BEC results from the bal- 
ance between the gravitational attraction and the quan- 
tum pressure arising from the Heisenberg uncertainty 
principle. Using dimensional analysis in Eq. (|27p . i.e. 

lrr?R^ ^ GM/R^, we obtain the length-scale 



Rq 



GMr 



(48) 



which gives the typical size of a self-gravitating BEC with 
mass A/ without short-range interaction (a = 0). 

In the TF approximation, in which the quantum po- 
tential is negligible, the structure of the self-gravitating 
BEC results from the balance between the gravitational 
attraction and the short-range repulsion due to scatter- 
ing (when a > 0). Using dimensional analysis in Eq. 
([27ll. i.e. {ah^ /m^R^){M/R^) ~ GM/R^, we obtain the 
length-scale 



Ra — 



ah'' 
G^ 



1/2 



(49) 



which gives the typical size of a self-gravitating BEC with 
scattering length a > in the TF approximation. 

Considering Eq. ([27]) again, the quantum pressure 
and the pressure arising from the short-range interac- 
tion become comparable when {ah^ /m^R^){M/R^) ~ 
n^/m^R'^, i.e. Na/R ~ 1. Estimating R by Eq. (gS]) 
or (I49p . this condition can be rewritten x ~ 1 where we 
have introduced the important dimensionless parameter 



(50) 



For a > and x ^ 1 , we are in the TF limit in which the 
quantum potential is negligible. In that case, the equi- 
librium state results from a balance between repulsive 
scattering and gravitational attraction. Alternatively, for 
X "C 1, we are in the non- interacting limit in which scat- 
tering is negligible. In that case, the equilibrium state 
results from a balance between quantum pressure and 
gravitational attraction. The transition between these 
two regimes occurs for x ^ 1- For a given value of the 
scattering length a, the TF limit is valid for M ^ Ma 
where 



VGr. 



(51) 



and the non-interacting limit is valid for M <C Ma- For 
a given value of the mass M, the TF limit is valid for 
aq where 



GAPm ' 

and the non-interacting limit is valid for a <C flg jl44l | 



(52) 



When a > 0, the gravitational attraction is necessary 
to balance the repulsive quantum potential and the re- 
pulsive short-range interaction. When a < 0, we could 
expect an equilibrium between the repulsive quantum 
potential and the attractive short-range attraction, in 
which gravitational effects are negligible. Using dimen- 
sional analysis in Eq. i.e. {\a\h^ /m^R'^){M/R^) ^ 
h^/m'^R'^, we obtain the length-scale 



K^N\a\, 



(53) 



which gives the typical radius of a non-gravitational BEC 
with attractive short-range interactions. However, as we 
shall sec, such equilibria arc unstable. 



H. Analogies between bosons and fermions 

The mean-field Gross-Pitaevskii equation ^ also de- 
scribes a gas of fermions when one takes into account the 
quantum potent ial Q arising from the Heisenberg uncer- 
tainty principle jl45| . This description goes beyond the 
TF approximation and can be useful to regularize the 
dynamics at small scales [i^. In the case of fermions, we 
must also take into account the quantum pressure arising 
from the Pauli exclusion principle. It can be calculated 
from the Fermi-Dirac distribution function at T = 0. 
In d dimensions, the equation of state is p ~ Kp^^'^/'^ 
where K = (l/(rf + 2))(rf/25'rf)^/'^(27rm Vm^+^/'^ for spin 
s = 1/2 fermions [ll7[. This is equivalent to a polytrope 
of index n = (1/2. This pressure term is the one that 
appears in the hydrodynamic equation (jl4p . Using Eq. 
p^ . it corresponds to an effective potential of the form 
h{p) = {d/2 + l)Kp^/'^ in the Gross-Pitaevskii equation 
(l6|). In particular, in d = 3 dimensions, the pressure is 
p = Kp^l^ with K = (l/5)(3/87r)2/3(27r;i)VmS/3 leading 
to a potential h = {h/2)Kp'^/^ and to a CP equation 



ifi 



at 2m 2 m 

(54) 

In comparison, a self-gravitating BEC with a potential 
([5]) is described by the CP equation 

ih^ ^ Aijj + m^Th + N M ^P- (55) 

ot 2m m 

Apart from the difference in the exponent, we see that the 
I ■01^ (quartic) self-interaction of bosons (p oc p^ cx IV'I"') 
plays a role similar to the Pauli exclusion principle for 
fermions, equivalent to a |?/j|-'^°/3 interaction (p cx p^^^ cx 
10/3 -J _ Indeed, a quartic self-interaction with a > 
is equivalent to an inter-particle repulsion that domi- 
nates over uncertainty pressure for N ^ 1. Similarly, 
in the case of fermions, the exclusion pressure dominates 
over uncertainty pressure for N ^ 1. In the Newtonian 
regime, and in the TF approximation, self-coupled boson 
stars arc equivalent to n = 1 polytropes and fermion stars 
(like classical white dwarf stars) to n = 3/2 polytropes. 
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The analogy between boson stars with a self-interaction 
and fermion stars takes even more sense in the relativistic 
regime [HI]. In that case, a self-coupled boson star has an 
equation of state p = e/3 = Kn^^^ (where e is the energy 
density and n is the particle density) correspond ing to a 
n = 3 polytrope like in the core of neutron stars [32h35I | . 

We note that the potential h[p) (x p^^"^ associated with 
fermions becomes equivalent to the potential h{p) cx p 
associated with self-coupled bosons when d = 2. In fact, 
the dimension d = 2 is a critical dimension [118| . When 
we consider a gas of repulsive (impenetrable) bosons, the 
potential h(p) = gp arising in the GP equation ceases to 
be valid for d < 2 (in d = 2 it remains marginally valid 
with logarithmic corrections). In particular, in d = 1, it 
is replaced by h{p) = tt^Tt? p^ /2m^ exactly like for spinless 
(s = 0) fermions. This is a mani festa tion of the boson- 
fcrmion duality in one dimension jll9l |. 



III. EXACT RESULTS 
A. The energy functional 

The total energy associated with the GPP system 
([7]) or, equivalently, with the quantum barotropic Euler- 
Poisson system (HU]), (|14p and ([7]) can be written 

Etot = Qc + Qq + U + W. (56) 

The first two terms correspond to the total kinetic energy 



e 



2m 



IVV-P dr. 



(57) 



Using the Madelung transformation, it can be decom- 
posed into the "classical" kinetic energy 



and the "quantum" kinetic energy 

1 



''Q = — I PQ dr. 

m 



(58) 



(59) 



Substituting Eq. (|T2|) in Eq. (|59|) . we obtain the equiva- 
lent expressions 



0Q = / VpAV^dr 



JiVV-prdr^^J^^dr. (60) 



2m2 

The third term is the internal energy 
U = i pj'?^dp,dv 



= j [ph{p)-p{p)] dr^ J Hip)dr, (61) 



where we have used Eq. (|15p and integrated by parts to 
obtain the second equality. For a polytropic equation of 
state (fT7|) . it takes the form 



U 



7-1 



p'^ dr 



1 



7-1 



pdr. 



(62) 



In particular, for the potential ([S]), using Eq. ([TB|) . we 

get 



U =^ 



2Trah 



p^ dr. 



(63) 



Finally, the fourth term is the gravitational energy 



W = ^ I p^dr. 



(64) 



The total energy per particle can be expressed in terms 
of the wave function as 



Etnt — 



where H 
written 



S"^ . 1 

IVV'P + H{Nm\if) + -TO$|i/)p 

2m 2 



dr. 



(65) 



H/N. Then, the GP equation dH) can be 



dt ~ sr 



(66) 



For a power-law potential h{p) — [7^7/(7- l)]p'^ , we 
have H{p) = [A7(7- l)]p^ so that H{Nm\ij\^) = 
with K = KN^-^rrif/{j - 1). 

It can be shown (see Appendix [E| that the total en- 
ergy Etot is conserved by the quantum barotropic Euler- 
Poisson system (or by the GPP system). The mass 
M = J pdr is also conserved. Therefore, a minimum of 
the energy functional Etot[Pi u] at fixed mass M is a non- 
linearly dynamically stable steady state of the quantum 
barotropic Euler-Poisson system (this follows from gen- 
eral results of dynamical stability jl20| ) . We are therefore 
led to considering the minimization problem 



min {E'tot [p, u] 



M}. 



(67) 



An extremum of energy at fixed mass is given by the vari- 
ational principle SEtot — (xSAl = where a is a Lagrange 
multiplier taking into account the mass constraint. Us- 
ing the results of Appendix |E1 this gives u = and the 
condition 



m<i> 



ih{p) 



2m y/p 



(68) 



This equation is equivalent to the steady state equation 
provided that we make the identification 



a — E/m. 



(69) 
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This shows that the Lagrange multipUer (chemical po- 
tential) in the constrained minimization problem (|68p is 
equal to the eigenenergy E by unit of mass (if we ap- 
ply the variational principle (p7|) at Eq. (p5)) we get Eq. 
((22)) ). On the other hand, considering the second order 
variations of energy, we find that the distribution is dy- 
namically stable iff 



8m2 



S^Etot = 



^'5 



dr - 



h'{p){Spfdr 



.3/2 



dp5^ dr 



8m2 



(Spf dr > 0, 



(70) 



for all perturbations that conserve mass: / Spdr ~ 0. 

Remark 1: If we plot a = E/m = dEtot/dM (conju- 
gate quantity) as a function of M (conserved quantity), 
we can determine the stability of the system by a di- 
rect application of the P oinca re theory of linear series of 
equilibria (see, e.g. HO, Il2l| for details). According to 
the Poincare theorem, a change of stability can only oc- 
cur at a turning point of mass or at a bifurcation point. 
Therefore, if we know a limit in which the configuration 
is stable, then we can use the Poincare theorem to de- 
duce the stability of the whole series of equilibria. We 
shall use this method in Paper II. 

Remark 2: In the TF approximation, the energy func- 
tional ([55| reduces to the standard Chandrasekhar en- 
ergy functional. It is well-known that a polytrope with 
index 7 > 4/3, including the polytrope 7 = 2 (n = 1) 
corresponding to Eq. p6p . is a minimum of energy at 
fixed mass. Therefore, it is nonlinearly dynamically sta- 
ble with respect to the barotropic Euler-Poisson system. 
Its linear dynainical stability can also be settled by using 
the Eddington |l22l| equ ation of pulsation or the Ledoux 
(123} criterion (see [ll3| and Appendix B of [ll]). 



B. The virial theorem 

From the quantum barotropic Euler-Poisson system 
([T0|. dH]) and d?]), we can derive the general time- 
dependent virial theorem (see Appendix |F|) : 

h = 2{Q, + QQ) + 3 j pdr + W, (71) 

where 



/ = I pr'^ dr, 



(72) 



is the moment of inertia. For a polytropic equation of 
state (fT7|) . we have the identity J pdr = (j — 1)U. More 
specifically, for the potential (O leading to Eq. ([T5|) . we 
get J pdr = U. In that case, the time-dependent virial 
theorem can be rewritten 



-/ = 2(6^ + Sq) + 3U + W. 



(73) 



At equilibrium (/ = Oc = 0), we obtain the time- 
independent virial theorem 

2eQ + 3U + W = 0. (74) 

On the other hand, the energy functional ([56| reduces to 

Etot ^Bq + U + W. (75) 

Finally, multiplying the steady state equation ([M)) by 
p and integrating over the configuration, we obtain the 
general identity 



ph dr + 2W ^ NE. 



(76) 



For a polytropic equation of state ([T7| . we find that 
J phdr = jU. More specifically, for the potential (O 
leading to Eq. ((T6)) . we get J phdr = 2U. In that case, 
Eq. (j76p can be rewritten 



eQ + 2U + 2W = NE. 



(77) 



In the non- interacting case (a — 0), the internal energy 
vanishes: U = 0. The three independent equations 
([75| and ^ reduce to 2eQ + W = 0, Etot = eQ + W 
and Qq + 2W = NE. From these equations, we obtain 
the relation 



Etot - InE. 



(78) 



This relation shows that the total energy Etot is not equal 
to NE, as we could naively believe |55|. The 1/3 factor 
was previously obtained by Membrado et al. |112l | from 
a different argument. We can check that the above rela- 
tions are satisfied by the different components (j30p - ([33)) 
of the energy. 

In the TF approximation, the quantum energy is ne- 
glected: 6q = 0. The three independent equations (|74l) . 
([75| and ^ reduce toW + W = Etot = U + W and 
2U + 2W = NE. From these equations, we obtain the 
relation 



Et. 



1 



-NE. 



(79) 



We can check that the above relations are satisfied by 
the different components (|35]), (gH), (gS]) and dM]) of the 
energy. 

Finally, in the non-gravitational limit, the potential 
energy is neglected: W = 0. The three independent 
equations ^ and jZll) reduce to 2eQ + 3U = 

0, Etot ^ Qq + U and Qq + 2U = NE. From these 
equations, we obtain the relation 



Et. 



-NE. 



(80) 



IV. THE GAUSSIAN ANSATZ 



To obtain the density profile of a self-gravitating BEG 
and the mass-radius relation, we need to solve the dif- 
ferential equation ([27]) expressing the condition of hy- 
drostatic equilibrium. This will be done numerically in 
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Paper II. However, it can also be useful to obtain ap- 
proximate analytical results. In that respect, we shall 
follow an approa ch si milar to the one devel oped by Stoner 
(1929,1930) [2IIT21. Nauenberg (1972) and Cha- 

vanis (2007) |26| in the case of classical and relativistic 
white dwarf stars. The idea is to prescribe an approx- 
imate density profile (characterized by its mass M and 
radius i?), compute the total energy Etot{R) and mini- 
mize it with respect to i?, for a given mass M, in order to 
obtain the equilibrium radius R{M). This method has 
provided very good approximations of the mass-radius 
relation of white dwarf stars. We shall see that it also 
provides a good approximation of the mass-radius rela- 
tion of self-gravitating BECs. In order to evaluate the 
quantum kinetic energy which involves density gradients, 
we shall make a Gaussian ansatz. The Gaussian ansatz 
is particularly accurate for small or moderate values of x 
for which the density profile extends to infinity (e.g., in 
the non- interacting case X = 0). By contrast, in the case 
X » 1, this ansatz is poor because the density profile 
approaches that of a ri = 1 polytrope which has a com- 
pact support (this is the exact solution in the TF limit 
X — >■ -|-oo). However, for the sake of illustration, we shall 
use the Gaussian ansatz for all configurations. 



A. The energy functional 

We shall calculate the energy functional ([55]) by mak- 
ing a Gaussian ansatz for the density profile 



p(r) = M 



3/2 



(81) 



The central density is p(0) = M/{'k^/^R^) and the corre- 
sponding rotation curve is 



Vcir) 



f GM 



1/2 r 



<0 



1/2 



(82) 



For comparison with the results of Paper II, it is conve- 
nient to introduce the radius containing 99% of the total 
mass. For the Gaussian density profile ([5T|) . we find that 
i?99 = 2.38167iT!. On the other hand, using Eq. (l5Tt . 
the quantum kinetic energy, the internal energy and the 
potential energy are given by 



Or 



Tt^M 



w 



GAP 

/ 

R 



with the coefficients 
3 



C 



(27r)3/2 ' 



2tt 



(83) 



(84) 



To compute the potential energy, we have used the for- 
mula W = —AttG J^°° p{r)M{r)r dr valid for a spher- 
ically symmetric distribution of matter, and integrated 



by parts. The moment of inertia is given by 

9 3 

I = aMR\ with a = -. 



(85) 



Finally, using the velocity profile (jASP of Appendix [Aj 
we find that the classical kinetic energy is given by 



1 fdR 
Gc = -aM — 

2 \ dt 



(86) 



From these expressions, the energy functional can 
be rewritten as a function of R and R (for a fixed mass 
M) as 



1 



Etr,t = -aM[ — 



dR 



n^M 2Trah^AP GAP 



c- 



m^R^ 



R 



(87) 



This can be interpreted as the total energy 
1 fdRV , , 



of a fictivc particle with mass aM and position R moving 
in a potential 



V{R) = cr 



h^AI 
im?R^ 



+ c 



2'Kah^AP GAP 



R 



(89) 



We shall come back on this mechanical analogy in Sec. 
HVCl The potential V{R) is plotted in Fig. [T]for different 
values of M (with a > and a < 0). These different cases 
are studied in the sequel. 
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FIG. 1: Potential V{R) for different values of M. The ra- 
dius and the mass are normalized by Ra and Ma (with a 
replaced by \a\) and the potential by Ea = GM^/Ra = 
n{GmY/'^ /\a\^'^ . Then, V{R) = aM/R^ ± 2'k(,M'^/R^ - 
vM"^ /R (with -I- for a > and - for a < 0). For a > 0, 
we have taken M = 1.2 and for a < we have taken 
M = 0.9 < M^ax and M = 1.2 > M^a^r where M,^ax = 1.085 
(see below). 
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B. The mass-radius relation 

A minimum of the energy functional Etot [p, u] , given 
by Eq. ([55]) . at fixed mass M determines a stable steady 
state of the quantum barotropic Euler-Poisson system 
(equivalent to the Gross-Pitaevskii-Poisson system). In 
the Gaussian approximation, we are led to determining 
the minimum of the function EtotiR, R), given by Eq. 

at fixed mass M. Clearly, we must have i? = so 
that a minimum of energy at fixed mass is a steady state. 
Then, we must determine the minimum of the potential 
energy V{R). Taking the derivative of V{R) with respect 
to R, we obtain 



V'iR) 



-2a 



3C- 



2Trah'^M^ GAP 



i?2 



(90) 



The critical points of the potential energy V{R), satisfy- 
ing dV/dR = 0, are solution of the equation 



M = 



2a 



GmTR 



V 1 



(91) 



The radius can be expressed as a function of the mass 
according to 



R=^- 



V GMm^ 



l±\l 



GttCi' GmM^a 



(92) 



with + when a > and ± when a < 0. This equation pro- 
vides an analytical approximate expression of the mass- 
radius relation of a self-gravitating BEG with short-range 
interactions. It can also be obtained from the equilibrium 
virial theorem (|74p by making the Gaussian ansatz (see 
Sec. HVCl) . 

Let us first consider asymptotic limits of this relation, 
(i) In the non-interacting case (a = 0), we obtain 



R 



2a 



GMr 



(93) 



-99 



The radius Rgg containing 99% of the mass is R 
8.955h^ /GMm^ . This can be compared with the exact 
resuh ([29]) giving E^f"'^* = 9.9h'^/GMm\ The agree- 
ment is fairly good, (ii) In the TF approximation (when 
a > 0), we get 



\ v J V Gm^ 



1/2 



(94) 



The radius is independent on mass. The radius 
i?99 containing 99% of the mass is given by Rgg = 
4.125(a?iVG'm3)i/2. This can be compared with the ex- 
act result dSS]) giving = 2.9m{ar? /Grn?Y/'^ . The 
agreement is less good than in the non-interacting case. 
The reason is related to the fact that the distribution 
psp has a compact support so that it is quite different 
from a Gaussian. The radius defined by Eq. rep- 
resents the minimum radius Rmin of the self-gravitating 



BEG for a given value of the scattering length a > 0. For 
R — > Rmim the mass tends to -l-oo. According to Eq. 
(j9T|) . it diverges like 



M 



1/2 



1 



\fGma 



1 



(95) 



(iii) In the non-gravitational limit (when a < 0), we get 
SttC M\a\ 



R 



(96) 



The radius Rgg, containing 99% of the mass is given by 
Rgg = 1.900Af |a|/TO. However, we shall see that these 
configurations arc always unstable. 

A critical point of ViR) is an energy minimum iff 
V" (K) > 0. Computing the second derivative of V{R), 
we obtain 



V"{R) = 6a 



m?R^ 



2TTah^AP GAP 
12C — 2iy- 



m^R^ 



i?3 



(97) 



Using the mass-radius relation (|9ip . the foregoing equa- 
tion can be rewritten 



R-^ \ vL 



yGm^R^ ) 



(98) 
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FIG. 2: M as a function of Rgg for given a > 0. The 
mass is normalized by Ada and the radius by Ra. Thus, 
M = 2aR/{uR^ - Q-kC,) with Rgg = 2.38167i?. The ra- 
dius is given as a function of the mass by i? = [a /vAi){l -\- 
-|- (3nC,uA'P /a'^). In the non- interacting hmit M — >■ 0, we 
get M ~ 2a /uR i.e. M ~ 8.955/i?99 and in the TF limit 
M +00, we get R -> {6-k(,/uY^'^ i.e. Rgg 4.125. All the 
configurations are stable. 



Let us first consider the case of repulsive short-range 
interactions (a > 0). The mass-radius relation is repre- 
sented in Fig. [21 There exists one, and only one, solution 
for each value of the mass and it is stable since, according 
to Eq. it is a minimum of energy {y"{R) > 0); see 

Fig. [T] The radius is a decreasing function of the mass. 
The non- interacting limit corresponds to M ^ Ala, i-e. 
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FIG. 3: M as a function of Rgg for given a < 0. The mass 
is normalized by Ma and the radius by Ra (where a is re- 
placed by \a\). Thus, M = 2aR/{vR^ + QtvC,) with ijgg = 
2.38167-R. The radius is given as a function of the mass by 
R = [(j/uM){l ± a/1 — 67r("!/Af2/(T2). In the non-interacting 
limit i? -foo, we get Af ~ 2a/vR i.e. M ~ 8.955/i?99 
and in the non-gravitational limit i? — >■ (unstable), we get 
M ~ o-7?/(37rC) i.e. ~ 0.5262ii99. There exists a max- 
imum mass Mmax ~ a I ^/WkCj) = 1.085 corresponding to a 
radius R» = {QiiC,/vY^'^ i.e. _R99 — 4.125. The configurations 
are stable for R> R, and unstable for R < R,. 



with the same mass. However, according to Eq. (|98p . 
only the solution with the largest radius R > R^, is sta- 
ble (local minimum of energy V"{R) > 0). The other 
solution is an unstable maximum of energy {V"{R) < 0); 
see Fig. [T] We can check that the change of stability 
{V"{R) = 0) occurs at the turning point of mass in the 
series of equilibria (A/'(i?) = 0), in agreement with the 
Poincare theorem. In the stable region, the radius is a de- 
creasing function of the mass. The non-interacting limit 
corresponds to Af ^ Ma ~ M^ax and R ^ Ra ^ R*- 
For M — > 0, the radius R — > +oo with the scaling (|93p . 
For M — >■ Mmax, the radius R tends to the minimum 
value i?*. The non-gravitational limit corresponds to 
Af — > and i? — > but these configurations are in- 
accessible sinc e they are dynamically unstable (energy 
maxima) jl46j |. If the system is initially placed on the 
unstable branch, it is expected to undergo gravitational 
collapse {R{t) 0) or to evaporate (i?(i) +oo), see 
Fig. [T] It may also relax towards the stable equilibrium 
state with a larger radius {R{t) — > i?s) provided that it 
is able to dissipate energy, e.g. by radiation. 



C. The virial theorem 



R Ra ^ Rmin • 

For Ai" 0, the radius R -J> +oo with 
the scaling The TF limit corresponds to > Af^, 
i.e. R ^ Ra ^ Rmin- For A/ +oo, the radius R tends 
to the minimum value Rmin given by Eq. (|94|) with the 
scaling (|95p. In the non- interacting case a = 0, the mass- 
radius relation is given by Eq. (|93p which determines a 
stable equilibrium state for any mass M . 

Let us now consider the case of attractive short-range 
interactions (a < 0). The mass-radius relation is repre- 
sented in Fig. 131 There exists a maximum mass 



1/2 



corresponding to the radius 



(99) 



\ V J 



Gm? 



(100) 



It may be noted that, contrary to the Chandrasekhar 
mass or to the Kaup mass, this maximum mass is a 
purely Newtonian result. For that reason, it is gen- 
erally very small (see Appendix IB 4[) . The approx- 
imate values Mraax = 1.085?i/-\/Gm|a| and i?gg = 
4.125(|a|?i^/G'm^)^/^ obtained with the Gaussian ansatz 
are in fairly good agreement with the exact re- 
sults M^^^^ = im2n/^Gm\a\ and {Rl^Y'""'^ = 
b.b{\a\h^ /Grrv^Y^'^ obtained numerically in Paper II. For 
M > Mmax , there is no solution (no critical point of en- 
ergy) and the system undergoes gravitational collapse to 
a black hole. For M < Mmax, there exists two solutions 



Using the Gaussian ansatz, the time-dependent virial 
theorem ([75)) can be written 



1 ..d^R^ nrfdR 
-aM—r- = aM — 



+3C 



dt 



2-Kah^M^ 
m^R^ 



2cr 



?Ym 

m^R^ 
GM^ 
R ■ 



Since 



fR^ 

dt^ 



= 2R- 



dt2 



(101) 



(102) 



we note the nice cancelation of terms in Eq. (|10ip leading 
to the final equation 



, „ h^M 



2'Kafi^M'^ GM^ 
3C TTTi 1^- 



i?2 



(103) 



The equilibrium virial theorem (d^R/dt^ ~ 0) returns 
the mass-radius relation (|9ip obtained from the condition 
dV/dR = 0. In fact, the time-dependent virial theorem 
(|103p can be written 



aM 



fR 
dt^ 



dV_ 
dR' 



(104) 



This equation describes the motion of a fictive parti- 
cle with mass aM and position i? in a potential V{R). 
Therefore, the total energy Etot = Qc + V given by Eq. 
is conserved 



dEt. 



dt 



dt 



(Oe + y) =0. 



(105) 
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In this mechanical analogy, a stable equilibrium state cor- 
responds to a minimum of V{R) as we have previously 
indicated. Alternatively, Eq. ()104|) can be viewed as the 
Hamilton equation of motion of the fictive particle asso- 
ciated with the Hamiltonian 



D. The pulsation equation 

To study the linear dynamical stability of a steady 
state of Eq. (|104p , we make a small perturbation around 
that state and write R{t) = R + €{t) where R is the equi- 
librium radius and e{t) ^ R is the perturbation. Using 
V'{R) — and keeping only terms that are linear in e, 
we obtain the equation 



0, 



where uj is the complex pulsation given by 

1 



aM 



V"{R). 



(106) 



(107) 



According to this equation, a steady state is linearly sta- 
ble iff > that is to say iff it is a (local) minimum of 
energy V{R). In that case, the system oscillates about 
its equilibrium value with a pulsation w. Otherwise, the 
perturbation grows exponentially rapidly with a growth 
rate A = -v/^ 



Using Eq. (|97| for V"{R) and compar- 



ing with Eq. 



we find that 



(108) 



In the non- interacting case [U = 0), using the virial the- 
orem ([74|) . we get = —W/I. In the TF approxima- 
tion (0Q = 0), using the virial theorem ([71)) . we obtain 
= —2W/I. This expression coincides with the Ledoux 
formula for a polytrope of index 7 = 2 [l23l |. In paper 
III, we show that this result can be generalized to an 
arbitrary polytropic index (see also Ref. [26j'). 

Using Eqs. (fTUT)) . ((98)) and (l9T|), the pulsation can be 
written in terms of the radius as 



2a 1 



(109) 



For a > 0, the pulsation-radius relation is plotted in 
Fig. m We note that > for any R so that the con- 
figurations are stable. When slightly perturbed about its 
equilibrium state, the system oscillates with a pulsation 
oj. For a < 0, the pulsation-radius relation is plotted in 
Fig. m We find that > for i? > R^ (stable) while 
< for i? < i?* (unstable). In the latter case, the 
perturbation grows with a growth rate A = \J —uP' . 

Let us derive asymptotic expressions of the pulsation, 
(i) In the non-interacting case, we obtain 



V GM 2a r 
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FIG. 4: uj^ as a function of Rgg for given a > 0. The pul- 
sation is normalized by uja = [GMa/ R^Y^'^ = Grr? jah and 
the radius by Ra. Thus, = (2cr/a_R*)(l %-kC,IvR^) j {\ - 
QivC,/uR^) with Rgg = 2.38167i?. In the non-interactine 
limit R — >■ +CXJ, we get uj'^ ~ 2cr/ai?^ i.e. uJ^ ~ 



32.18/i?gg and in the TF limit 

{aV^/l9.a-K^(')l{R/Rrmr^ " 1) ^ 
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FIG. 5: Lj^ as a function of 7?gg for given a < 0. The pul- 
sation is normalized by uja = {GMa/ RtY^'^ = Gm'^ /\a\h 
and the radius by Ra (with a replaced by |aj). Thus, uJ^ = 
(2cr/ai?*)(l-67rC/i^i?^)/(l + 67rC/!^i?^) with i?gg = 2.38167i?. 
In the non-interacting limit R — >■ -|-oo, we get ui^ ~ la jaR^ 
I.e. Lo ~ 32.18/ii|9 and in the non-gravitational limit i? — >■ 
(unstable), we get uP' ~ —lajaR^ i.e. cj^ ~ — 32.18/i?99. 
The pulsation is maximum for _Rgg = 5.247 with the value 
= 10.02 10"^ 



When a 7^ 0, this expression is asymptotically valid for 
R — > +0O and M — > 0. (ii) In the TF approximation (for 
a > 0), we get 



2vGM _ _^v^^G^^^Mnf^ 
aR^ ~ a(67rC)3/2 a^/^n^ 



(111) 



This expression is valid for R — > Rmim he. M — > +oo, 
with M[R) given by Eq. (|95|) . In that limit w — > +oo. 
(iii) In the non-gravitational limit (for a < 0), we find 
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FIG. 6: Collapse of the BEC for M > M^ax in the case of an 
attractive self-interaction (a < 0). We have represented the 
BEC radius R{t) as a function of time, by solving Eq. (|103p 
obtained with the Gaussian ansatz. The radius is normalized 
by Ra, the mass by Ma and the time by ta ~ \a\h/Gm'^ . The 
equation of motion can then be rewritten aMR — —V'{R) 
with ViR) = aM/R^ - 2TrCM^/R^ - vM^/R. The solu- 
tion is (2/aM)^/2f = dR/^V{Ro) - V{R) (solid line). 
The collapse leads to a finite time singularity. The col- 
lapse time tcou is obtained from the foregoing expression 
by using R{tcoii) = 0. The solution can then be rewritten 
{2/aMy^\t,ou -t) = /g"'*' dR/y/ViRo)-V{R). For t ^ 
tcoii, the radius behaves like R{t) ~ {25TrCM/ay/-' {tcou-tf^^^ 
(dashed line). We have taken A4 = 1.2, Rq — 1 yielding 
{2/aAdy/h,oii = 1.692. 



that 



2a fi^ 



(112) 



This expression is valid for i? — > and A/ — > 0. 

Remark 1: in the TF limit (when a > 0), the density 
profile is known analytically, see Eq. (pS]). Ther efore, we 
can be more precise. If we use the Ledoux |123{ formula 
up = —2W/I for the pulsation (which is approximate) 
with the exact expressions (|15|) and (|^T1) of W and /, we 
obtain u^Ledoux = 0.3512{GM/Riy^^. The Gaussian ap- 
proximation yields ujoauss = 0.3199{GM/R^y/^ corre- 
sponding to Eq. (|llip . Fi nally , the exact result obtained 
by solving the Eddington |l22| equation of pulsation nu- 
merically is ujexact = 0.3480(GM/i?3) . We note that 
the reference pulsation can be written {GM / R^Y^'^ ~ 

^l/2g5/4^11/4^-3/4^-3/2_ 

Remark 2: in the case of an attractive self-interaction 
(a < 0), there is no equilibrium when M > Mmax and 
the system is expected to collapse. We could attempt 
to describe this collapse by using Eq. (|103p as done 
in Fig. ini For R{t) 0, it reduces to aMdpR/dt^ ^ 
Q-nQah^M"^ /rrt'R'^ so that gravitational effects become 
negligible. This leads to a finite time collapse with a 
radius scaling like R{t) cx {tcoii ~ t)^/^. However, this 
scaling is different from the scaling R{t) cx {tcoii — t)^^"^ 



obtained by directly s olvin g the Gross-Pitaevskii equa- 
tion (without gravity) jl26l |. This shows that the Gaus- 
sian ansatz may not always be accurate. Eq. (|103p can 
be used close to a steady state so as to provide a good 
approximation of the pulsation period and of the growth 
rate, but it may lead to inaccurate results in more general 
situations. 



E. The radius versus scattering length relation 

For a given mass M, Eq. determines the radius 

i? as a function of the scattering length a. Inversely, Eq. 
(PT|) yields 



^ , vGMR-' -2<7-!—R 



(113) 



On the other hand, using Eqs. (fT07| . (|98|) and (fm]) . the 
pulsation is given in terms of the radius by 



2vGM 
ai?3 



1 - 



'GMm^R 



(114) 



The radius versus scattering length relation is repre- 
sented in Fig. [7] and the pulsation-radius relation is rep- 
resented in Fig. [5] 

Let us first consider the case of repulsive short-range 
interactions (a > 0). There exists one, and only one, so- 
lution for any value of a and it is stable (w^ > 0). In 
the non- interacting case a ~ 0, the radius is given by 
Eq. dM]) and the pulsation by Eq. pTO)) . In the TF 
approximation, valid for a 3> ag, the radius is given by 
Eq. (|94l) and the pulsation by Eq. (jllip . The radius in- 
creases with the scattering length a and tends to -|-oo for 
a — > +00. Therefore, a repulsive self-interaction allows 
to construct dark matter halos whose size is much larger 
than for systems without self-interaction. 

Let us now consider the case of attractive short-range 
interactions (a < 0). There exists a minimum scattering 
length 



corresponding to the radius 



i?* — — 



1/ GMm^ 



(115) 



(116) 



The approximate values an 



-l.USh^/GAPm and 



i?gg = 4.477?i /GMm obtained with the Gaussian 
ansatz are in fairly good agreement with the exact 
results a^™'^* = -imdh^ /GAPm and (RggT'""'^ = 
b.Qh^ / GMm? obtained numerically in Paper II. For a < 
CLmim there is no equilibrium and the system undergoes 
gravitational collapse. For a„j„ < a < 0, there exists 
two solutions with different radii. The solution with the 
largest radius i? > i?* is stable (w^ > 0) and the solution 
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with the smahest radius i? < i?* is unstable (w^ < 0). In 
the stable region, the radius increases with the scatter- 
ing length a. There exists a minimum radius i?* given 
by Eq. (|116p . For i? — > i?*, the pulsation behaves like 
uj^ = {2vGM/aRt){R - i?*). In the non-gravitational 
limit, corresponding to i? ^ i?*, the radius is given by 
Eq. dnni) and the pulsation by Eq. ([TT^ . 
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FIG. 7: Radius Rgg versus scattering length a relation for a 
given value of the mass M. The radius is normalized by Rq 
and the scattering length by ag. Thus, a — (vR? — 2gR) /67r(" 
with Tigg = 2.38167-R. The radius is given as a function of the 
scattering length hy R = {g / a/1 + &TiC,i'a/ a^) with the 
signs -I- and — for amin < a < and with the sign -I- for a > 0. 
Stable solutions exist for a > amin = —a'^ /{6it(^u) = —1.178 
and Rgg > i?gg corresponding to R, = a/v i.e. i?gg = 4.477. 
The radius R is monotonically increasing with the scattering 
length a. The non-interacting limit (a = 0) corresponds to 
R = 2cr/i/, i.e. 7?gg = 8.955. In the TF limit, valid for 
a +00, we have R ~ {6TvC,/uy^^a^^^ , i.e. i?gg ~ 4.1252a^/^ 
The non-gravitational limit, valid for R (lower branch) 
corresponds to i? ~ {3Tv(^/a)\a\, i.e. -Rgg ~ 1.900|a|, but these 
solutions are unstable. 
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FIG. 8: Pulsation oj^ versus radius i?gg relation for a given 
value of the mass M. The pulsation is normalized by 
UQ = (GM/R^y^'^ = G'^M^m'^/h^ and the radius is nor- 
malized by Rq. Thus, tj^ = {2u/aR^){l - a/uR) with 
i?gg = 2. 3816771. In the non- interacting case (a = 0), 
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In the TF limit R -i> -|-oo, 
uj- ~ 2u/aR-^ ~ 7.186/i?gg. In the non-gravitational limit 
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-32.18/7?g9. The pulsation IS max- 
8.443 10"^ The 



-R — >■ 0, a; ~ — 2(7/q_R ~ — o^.±o/-n,gg. 
imum for i?gg = 5.970 with the value 
pulsation can be expressed in terms of the scattering length 
as = ±{2v'^/aa^)^l + Q-KC^va/a'^ /{I ± -s/l + &Ti(,va/ o'^Y . 



infinite system (the case of an expanding universe is 
treated inJlM), it is necessary to modify the Poisson 
equation [l28l |. This can be done in different ways, either 
by subtracting the average density p to the local den- 
sity p{r,t) or by introducing a shielding length k,~^ in 
the interaction and letting k 0. This is not a "Jeans 
swindle" as oftentimes said |llC| but rather a well-defined 
and rigorous mathematical procedure ^128j . We therefore 
consider the system of equations 



|+V.(pu)=0, 



(117) 



V. JEANS-TYPE INSTABILITY OF A 
SELF-GRAVITATING BEG 

In this section, we study the linear dynamical stabil- 
ity of an infinite homogeneous system of self-gravitating 
BECs described by the quantum barotropic Euler equa- 
tions. This is a generalization of the classical Jeans prob- 
lem to the fully quantum context. This type of anal ysis 
has been , performed by several authors [7l|, Izl: US lOOl - 
Il02l . Il27l | in different contexts. However, these authors 
did not take into account short-range interactions that 
play an important role in the physics of cosmic BECs. 
Therefore, our discussion is more general. We shall also 
develop the connection with the results established in the 
previous sections. 

Let us consider the quantum barotropic Euler-Poisson 
system (|10p . ((T?)) and ([7]) describing the dynamical evo- 
lution of an infinite system of self-gravitating BECs. In 
order to correctly define the gravitational force in a static 



du 

at 



+ (u-V)u 



1. 



(119) 



where p is the average density. In our study, the temper- 
ature r = and the pressure p{r, t) is due to short-range 
interactions (scattering) between particles. For specific 
applications, we shall assume a barotropic equation of 
state of the form p^ . 

The linearized quantum barotropic Euler equations 
around an infinite homogeneous distribution with p = p, 
u = and $ = are 



(120) 



dSu 



Am? 



V{ASp), (121) 
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A5$ = inGSp, (122) 

where = p'{p) = ph'{p) is the velocity of sound. For 
the equation of state ([T6l) . it is given by 



Anah p 



(123) 



It is easy to combine these equations into a single equa- 
tion for the perturbed density 5p. We find that 



d'^5p 



Am? 



A'^6p + clMp + -i-nGpSp. 



(124) 



Expanding the perturbation in plane waves of the form 
6p{r,t) oc cxp[i(k • r — cot)], we obtain the dispersion 
relation 



4^2 



c^k^ - AnGp. 



(125) 



This is the gravitational analogue of the Bogoliubov [12£ 
energy spectrum of the excitation of a weakly inter- 
acting Bose-Einstein condensate. For large wavenum- 
bers (small wavelengths) , the quasi-particle energy tends 
to the kinetic energy of an individual gas particle and 
w ~ hk? /2ra. The group velocity is 



dk 



2m2 



:fc. 



4m^ 



c2fc2 - AnGp 



(126) 



In the non- interacting case (a = 0), the pressure is zero 
(j) = 0) and the particles interact only via gravity. The 
dispersion relation reduces to 



Am? 



AirGp. 



(127) 



This equation exhibits a characteristic wavcnumber due 
to the interplay between gravity and quantum effects 
(Heisenberg's uncertainty principle): 



IGnGm'^i 



1/4 



(128) 



This quantum Je ans scale appears in the works of [71|, |72|, 
IsqI . IqoI . I102| . Il27l | . The system is stable for perturbations 
with k > kj and unstable for perturbation with k < kj. 
The maximum growth rate corresponds to fc = leading 
to 7 = \/ A'KGp. The mass contained within the sphere of 
diameter Aj, where Aj = iir/kj, is 



Ml 



3/4 



■,1/4 



(129) 



We therefore expect that the gravitational collapse of 
a homogeneous distribution of noninteracting bosons at 
T — Q leads to objects with typical radius Rj = Xj/2 
and typical mass Mj , or larger (recall that the maxi- 
mum growth rate corresponds to A — > -l-oo) |147| . The 



physical mechanism that leads to a non- vanishing Jeans 
scale and Jeans mass has the same nature as that which 
accounts for the equilibrium of the boson stars studied 
by Rufhni & Bonazzola [5^1 • It corresponds to a balance 
between the gravitational force and the quantum pres- 
sure. Eliminating the density between Eqs. p28p and 
(|129p . we obtain the formula 



12 Gm^Rj ' 



(130) 



which qualitatively agrees with the mass-radius relation 
obtained in [ssj . 
In the Thomas-Fermi approximation in which the 
quantum potential can be neglected, the particles inter- 
act via gravity and they experience a pressure due to 
short-range interactions. The dispersion relation reduces 
to 



c2fc2 - AnGp. 



(131) 



This is the usual Jeans dispersion relation. For a < 0, 
the system is always unstable. For a > 0, the Jeans 
wavcnumber is 



k. 



^/4^ fGm^ 



1/2 



(132) 



The system is stable for perturbations with k > kj and 
unstable for perturbation with k < kj. The maximum 
growth rate corresponds to fc = leading to 7 = y/AnGp. 
The characteristic wavcnumber (jl32p arises due to the in- 
terplay between gravity and scattering. We note that the 
Jeans wavcnumber is independent on the density. The 
mass contained within the sphere of diameter Aj, where 
A J — 27r/fc,7, is 



3/2 



TT / Aii'^ah? 
6^ I Gm3 



3/2 



(133) 



We therefore expect to form objects with typical radius 
Rj = Aj/2 and typical mass Mj, or larger. The physical 
mechanism that leads to a nonvanishing Jeans scale and 
Jeans mass has the same nature as that which accounts 
for the equilibrium of self-gravitating BECs with repul- 
sive short-range interactions studied by Bohmer & Harko 
[9^. It corresponds to a balance between the gravita- 
tional force and the pressure due to repulsive scattering. 
In fact, Eq. (|132p agrees with the radius of boson 
stars in the TF limit and Eq. (|133p corresponds to the 
mass-central density relation (|40p . 

In the non-gravitational case, the dispersion relation 
reduces to 



2 7,4 



rfc 
Airi? 



(134) 



For a > 0, the system is always stable. For a < 0, the 
critical wavcnumber is 



167r|a|p 



1/2 



(135) 
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The system is stable for perturbations with k > kj and 
unstable for perturbation with k < kj (the subscript J 
- for Jeans - is here an abuse of language since gravity 
is neglected). This characteristic wavenumber arises due 
to the interplay between quantum pressure and attrac- 
tive scattering. The mass contained within the sphere of 
diameter A,/, where Xj — 27r/fc,/, is 



Mt = - 



1 / 7rm 



3/2 



(136) 



We therefore expect to form objects with typical radius 
Rj = A,//2 and typical mass AIj, or larger. Eliminating 
the density between Eqs. (|135p and (|136p we obtain 



Mj = 



12|a| 



Ri 



(137) 



which reproduces the scaling of Eqs. ([55]) and 
The growth rate is maximum for k^ — (87r|a|p/7Ti)^/^ and 
its value is 7* = 47r|a|7ip/?7i^. 
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FIG. 9: Dimensionless dispersion relation fi^ = + 2aK^ — 
1 with SI — Lu/oJo, K = k/ko and a — a/ao where ljo ~ 
^/4^Gp, fco = (le-rvGpm^/h'^y^'^ and gp = {G m^ /Aiyph^y^^ . 
The Jeans wavenumber is Kj = —a + y/a'^ + 1. For a < 0, the 
maximum growth rate is 7, — \/l + a'^ reached for k» = —a. 
The TF limit corresponds to -ti the non-gravitational 
limit to 3> and the non-interacting limit to \a\ <c; 1. 



From the dispersion relation (|125|) . we note that the 
TF approximation is valid for k ka with ka = 
(167r|a|/9/7Ti)^/^ and the non-gravitational approximation 
is valid for k:^ kg with kg = {Gw? /\a\fi^Y/'^ . The non- 
interacting limit corresponds to ka ^ kg i.e. a ^ oq = 
now consider the general case. 
The pulsation vanishes at the critical Jeans wavenumber 



2m-' 




(138) 



Using the expression (|123p of the velocity of sound, this 



can be rewritten 



87r|a|p 



Awa^h p 



sgn(a) 



(139) 



The system is stable for k > kj and unstable for k <kj. 
From the Jeans wavenumber (|139p . we can define the 
Jeans radius Rj = Aj/2 = 7r/fcj and the Jeans mass 



Mj = |7ri?}. We find that 



Ml 



P 



1/2 



1 



+ 4^ - Sg-(«) 



3/2 ■ 



(140) 



Eliminating the density between Eqs. (|139p and (jl40p . 
we obtain 



Mj 



tRj 



12 GRj 



(141) 



This expression returns the mass-radius relation (j9ip . 
For a > 0, u!^ increases monotonically with k. Accord- 
ingly, the growth rate is maximum for k ~ (infinite 
wavelength) leading to 7 = ^/AnGp. For a < 0, 
achieves a minimum negative value at 



k., = 



87r|a|/9 



1/2 



(142) 



Accordingly, the growth rate is maximum for k = fc* and 
its value is 



> = V ^4 +4-Gp. 



(143) 



For 



7* 



> {G m^/Air h'^p)^/^, we find that 
ATrlaslTip/m^ ^ y/AnGp. Therefore, an attractive short- 
range interaction (a < 0) increases the growth rate of the 
Jeans instability. The pulsation is plotted as a function 
of the wavenumber in Fig. ^ for positive and negative 
scattering lengths. 

These results may have profound cosmological impli- 
cations. Dark matter is usually described by hydrody- 
namical equations of the form (|117p - (|119p without the 
quantum potential ((5 = 0). In the context of cold dark 
matter (CDM) models with vanishing temperature and 
pressure (T = p = 0), the usual Jeans analysis predicts 
that all scales are unstable. Consequently, the Jeans 
scale and the Jeans mass are zero (A,/ = Mj = 0). 
This is the intrinsic reason why CDM models generate 
cuspy dark matter halo profiles and an abundance of low 
mass halos. However, these cusps and satellites are not 
seen in observations [tJ [t^ ■ As argued by several au- 
thors t9fl, .102„ . ,127} , these problems could be solved if the 
dark matter of the universe is a self-gravitating BEC. In 
that case, the wave properties of the dark matter can 
stabilize gravitational collapse, providing halo cores and 
suppressing small-scale structures. Indeed, if dark mat- 
ter is a BEC, even at T = there exists a non-zero 
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Jeans length (Aj > 0). Stability below the Jeans scale 
is guaranteed by the Heisenberg uncertainty principle. 
Therefore, there exists a minimum radius and a minimum 
mass at which the system is stable. If the particles have 
a repulsive self-coupling (a > 0), stability results from 
the pressure arising from the repulsive scattering j95l |. 
This non-thermal pressure indeed stabilizes the system 
against gravitational collapse and leads to dark matter 
halos with a flat core equivalent to n = 1 polytropes 
(or other barotropic structures depending on the form of 
the self- interaction). Alternatively, if the particles have 
an attractive self-coupling (a < 0), generating a nega- 
tive pressure, this (anti)-pressure can enhance the Jeans 
instability and fasten the formation of structures. The 
virtues of these results could be combined by assuming 
that the sign of the scattering length a changes in the 
course of the evolution. It could be negative in the early 
universe to help with the formation of structures and 
become positive once the structures start to form in or- 
der to prevent their complete gravitational collapse. The 
mechanism of this change of sign is, however, unknown so 
that this idea remains highly speculative. In any case, the 
quantum pressure and the scattering pressure are small- 
scale effects. The BEC model and the CDM model differ 
at small scales only while they are indistinguishable at 
large scales [o^l- Therefore, at large scales of observa- 
tional interest, we recover the usual equations of CDM 
which have proven to be very relevant. 



VI. CONCLUSION 

In this paper, we have studied the structure and the 
stability of a self-gravitating BEC with short-range in- 
teractions. We have connected the results of Rufhni & 
Bonazzola [s^ obtained in the absence of self-coupling 
to the results of Bohmer & Harko [9^1 obtained for self- 
coupled BECs in the Thomas-Fermi approximation. We 
have also considered the case of attractive short-range in- 
teractions (a < 0) and found the existence of a maximum 
mass above which no equilibrium state exists. These re- 
sults have been obtained analytically by using a Gaussian 
ansatz and developing an analogy with a simple mechan- 
ical problem. In Paper II, we shall compare our approxi- 
mate analytical results with the exact results obtained by 
numerically solving the equation of hydrostatic equilib- 
rium. Finally, in Paper III, we shall extend our analytical 
method to more general situations. 

Our study was motivated by the proposal that dark 
matter halos could be gigantic cosmic BECs [tI, [13, l90l - 
[95| . In that case, gravitational collapse is prevented by 
the Heisenberg uncertainty principle or by the short- 
range interaction. This suggestion remains, of course, 
highly speculative since we do not know the true nature 
of dark matter. There exists other theories according to 
which dark matter could be made of massive neutrinos 
[Sjsoj . In such theories, gravitational collapse is pre- 
vented by the Pauli exclusion principle for fermions. On 



the other hand, whatever the nature of its constituents, if 
we view dark matter as a coUisionless system described by 
the Vlasov equation, dark matter halos could result from 
a process of violent coUisionless relaxation [s^ [s^ . In 
that case, gravitational collapse isprevented by Lynden- 
BcU's type of exclusion principle |5l| . This form of relax- 
ation is much more rapid and efficient than a "collisional" 
relaxation (e.g. for fermions). Furthermore, it generates 
a density profile with a flat core (containing possibly a 
massive degenerate nucleus) and a density halo yield- 
ing flat rotation curves. These features are remarkably 
consistent with observations making this alternative sce- 
nario quite attractive. It may also be recalled that the 
very existence of dark matter is questioned by some au- 
thors who intro duc ed modified gravity theories, like the 
MOND theory (130| |, to explain the astrophysical obser- 
vations without invoking dark matter. 



Appendix A: The velocity field 

In this Appendix, we determine the expression of the 
velocity field used in Sec. IIV Al to compute the classical 
kinetic energy ((55)) . 

The continuity equation (jlQ[) can be rewritten 



i91np 

dt 



V-u-|-Vlnp-u = 0. 



(Al) 



Let us assume that /9(r, t) is given by Eq. ([8T|) where 
R = R{t) is a function of time. Then, we have 



dlnp dlnR d\nR 

= -3— — + 2- 



dt 



dt 



i?2 dt 



(A2) 



(A3) 



Assuming that the velocity profile is of the form u(r, t) ^ 
f{t)r and substituting these relations in Eq. (jAl[) . we 
obtain 



dlnR 
dt 



fit) 



0. 



This relation is identically satisfied if f{t) 
Therefore, we can make the ansatz 



u(r,t) 

for the velocity field. 



dhiR 



dt 



R 

-' = r'^ 



(A4) 
dlnR/dt. 

(A5) 



Appendix B: Relativistic extensions and maximum 
mass 

In this Appendix, we propose a simple way to general- 
ize our Newtonian results to the relativistic regime. Our 
approach returns the right scaling of the maximum mass 
in the known cases and provides interpolation formulae 
in more general cases. 
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1. Thomas-Fermi limit for fermion stars 

The mass-radius relation of a non-relativistic fermion 
star is 



(Bl) 



with a prefactor 4.51 According to this relation, there 
exists an equilibrium state for any value of the mass M. 
In reality, this relation ceases to be valid when the radius 
approaches the Schwarzschild radius Rs = 2GM /c^ so 
that general relativistic effects come into play. Equating 
the two relationships, and introducing the Planck mass 
Mp = (hc/Gy^^, we obtain the maximum mass of a 
relativistic fermion star 



M, 



Ch 



Ml 



(B2) 



Therefore, the radius of a relativistic self-gravitating 
EEC without self-interaction scales like the Compton 
wavelength of the particles that compose the EEC. More 
precisely, the radius containing 95% of the mass is given 
by Eq. with a prefactor 6.03 [6l|. We note that the 
ratio between the classical radius (|E4p and the Compton 
wavelength scales like 



R 



Kaup 



M 



(E7) 



with a prefactor 15.6. Therefore, for non-relativistic 
EECs, such as dark matter galactic halos, the typical 
radius of the halo is much larger than the Compton wave- 
length (i? > Ac) since M < Mxaup- 



Thomas-Fermi Hmit for boson stars 



This corresponds to the Chandrasekhar mass which 
scales like Eq. (jE2[) with a prefactor 0.376 in general rel- 
ativity [3^ . The corresponding radius, which can be in- 
terpreted as the minimum radius of a relativistic fermion 
star is 



Mp 



Xr 



(B3) 



where Ac = fi/mc is the Coinpton wavelength. The pref- 
actor in Eq. (|B3]) is 3.52 [s^- Since m/Mp < 1, we 
note that i?mm ^ Ac. We also note that Ac/^_p — Mp/m 
where Ip = (hG/c^)^/'^ is the Planck length. 



2. Non-interacting boson stars 

In the absence of short-range interaction, the mass- 
radius relation of a non-relativistic self-gravitating EEC 
is given by [ssllTlll: 



R' 



GMm? ' 



(E4) 



with a prefactor 9.9 (if R represents the radius contain- 
ing 99% of the mass). This relation is valid as long as 
the radius is much larger than the Schwarzschild radius 
Rs = 2GM/c^. Equating the two relationships, and 
introducing the Planck mass, we obtain the maximum 
mass of a relativistic self-gravitating EEC without self- 
interaction 



M 



Kaup 



Ml 



(E5) 



This corresponds to the Kaup [SJ] mass which scales like 
Eq. (jE5[) with a prefactor 0.633. The corresponding 
radius, which can be interpreted as the minimum radius 
of a relativistic self-gravitating EEC without short-range 
interaction is 



R-n 



Ac. 



(E6) 



Wc now consider self-gravitating EECs with repulsive 
self- interaction (a > 0). In the TF limit, the radius of 
the configuration is given by 



Gm- 



1/2 



(E8) 



with a prefactor tt [95|- It is independent on the total 
mass M of the system. In order to make contact with 
studies that consider self-coupled particles interacting via 
a \X\(j)\^ potential [sl [H |9J, lill , we introduce the di- 
mensionless parameter 



Stt 



a 



(E9) 



measuring the strength of the short-range interaction. 
For a = 10^ fm and m = 1.44 cV/c^ wc get A/Stt = 0.01 
and for a = 1 fm and m = 14meV/c^ we find that 
A/Stt = 10^^°. In terms of this parameter, the radius 
can be rewritten 



R 



/ xh^y/' 1 

\8ttGcJ m 



(BIO) 



with a prefactor tt. Interestingly, the same relation ap- 
pears in the work of Arbcy et al. [93 | who consider a 
self-coupled charged scalar field and obtain an equation 
of state p = {Xfi^ /Am^c)p'^ equivalent to Eq. ([TH). The 
same equation of state is obtained by Colpi et al. [66j 
at low densities while p — p(? at high densities (see 
Appendix |D|) . There exists therefore a close connection 
between a self-coupled EEC described by the CP equa- 
tion and a self-coupled charged scalar field. Using Eq. 
(|E9|) . the parameter x introduced in Sec. IIIGI can be 
written x = (A/Stt) (M/Mp)^. The TF approximation is 
vahd provided that (A/Stt) (M/Mp)^ > 1. 

The Newtonian approximation is valid as long as the 
radius (|ESI) is much larger than the Schwarzschild radius 
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Rs = 2GM /c?. Equating these two relationships, we 
obtain the maximum mass 



M 



(Gm)3/2 



(Bll) 



Using Eq. ()B9p and introducing the Planck mass, it can 
be rewritten 



M 



(B12) 



This is the maximum mass appearing in the work of Colpi 
et al. with a prefactor 0.062. As emphasized by 

these authors, for A ~ 1, it scales like the Chandrasekhar 
mass Mp/vn? while the Kaup mass scales like Mp/m. In 
the presence of self-interaction, the maximum mass of a 
relativistic BEC is much larger than the Kaup mass by 
a factor Mp/m ^ 1, so that it becomes astrophysically 
relevant. The radius of the configuration is still given, in 
the relativistic regime, by an equation of the form (jB8[) 
or (jB10[) but the prefactor is different. It can be written 



m 



(B13) 



with a prefactor 0.3836 [131| . Therefore, for A ~ 1, the 
radius of a self-coupled BEC is much larger than the 
Compton wavelength since Mp/m ^ 1. 

More generally, the mass (|B12p of a self-coupled BEC 
is much larger than the Kaup mass (jBSP provided that 
A ^ (m/Mp)^. This is easily realized, even for weak 
self-interactions, since m/Mp <^ 1. Following Colpi et 
al. [6^, it makes sense to introduce the parameter (we 
have not written their factor 47r): 



A 



A^. 



Then, Eqs. (|B12p and (|B13|) can be rewritten 



M 



R 



(B14) 



(B15) 



with prefactors 0.062 and 0.3836. The self-interaction 
is important for A ^ 1 while it is negligible for A <C 1. 
Even for small A, the self- interaction is important because 
A can be quite large due to the greatness of the term 
(Afp/m)2. 

Finally, we emphasize that the radius i? of a self- 
coupled Newtonian boson star (see Eqs. (|B8|) and (|B10|) ) 
depends on m and a (or A) only through the combination 
a/m^ (or A/m^). The same observation holds for the ra- 
dius R and the mass M of a self-coupled relativistic boson 
star (see Eqs. (|BTT|) . ((BT21) and ((BTS)) ). 



4. Boson stars with attractive self-coupling 

For self-gravitating BECs with attractive self-coupling 
(A < 0), using Eq. (|B9[) . the maximum mass ([M)) and 



the minimum radius (jlOOp can be written 



(B16) 



with prefactors 5.073 and 1.1 (if R represents the radius 
containing 99% of the mass). For |A| ~ 1, M,nax is of 
the order of the Planck mass Mp = 2.1810~*kg, i.e. 
ridiculously small. This essentially means that a self- 
gravitating boson gas with attractive interactions is very 
unstable. However, the mass increases when |A| — )■ 
while the radius decreases. The system becomes rela- 
tivistic when the radius R approaches the Schwarzschild 
radius Rs = 2GM/c^ defined with the mass M. This 
happens when |A| < {m/Mp)"^ yielding a maximum mass 
M ^ Mp/m, equivalent to the Kaup mass, and a radius 
R ^ Ac of the order of the Compton wavelength. There- 
fore, when self-coupling is attractive, the maximum mass 
is very low (implying a strong instability) except if the 
mass of the bosons, or their self-interaction, is extremely 
small (which is possible for axions and ultra-light bosons 
ilQ, 65, 90J). 

Summarizing, for A > (m/Mp)^ the maximum mass 
scales like \f\M\/m^ [6^1, for A = it scales like M|,/m 
[Hi] and for A < ~{m/Mpf it scales like Mp/VX. It 
is interesting to obtain the three mass scales Mp/m^, 
Mp/m, Mp in the different regimes of a self-gravitating 
BEC. 



5. General case: an interpolation formula 

In the general case, the mass-radius relation of a non- 
relativistic self-gravitating BEC with short-range interac- 
tions can be approximated by the relation (see Eq. (|9ip ): 



M 



Gm^R 



1 - 



(B17) 



where we have get rid of (uncertain) numerical factors. 
In the relativistic regime, equating the radius R with the 
Schwarzschild radius Rs — 2GM /(? and using Eqs. (jB9p 
and (|B14p . we obtain after simplification 



M 



A, 



R - Ac\/1 + A. 



(B18) 



These formulae are valid for A > —1 (i.e. A > 
— {m/MpY). For A < — 1, the maximum mass and the 
minimum radius are given by the Newtonian expressions 
(|B16p . In the non- interacting case A = and in the TF 
limit A — > +00, we recover the previous scalings (|B5p . 
(jB6p and (jBlSp . Taking into acco unt the correct pref- 
actors obtained in [13, [6lj, [H, Il3l| , we can propose the 
interpolating formulae 



M = 0.633A/r 



A#2 

0.0096A— ^, 
m 



(B19) 
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i?95 = e.OSVl + 0.004AAc 



(B20) 



for the maximum mass and the minimum radius of a 
relativistic self-gravitating BEC with short-range inter- 
actions. Interestingly, an expression for the critical mass 
of a relativistic boson star with repulsive self-interaction 
similar to Eq. (jB19P has been obtained in [l^l based on 
different arguments. 



Appendix C: Jeans instability with a cosmological 
constant 

In this Appendix, we extend the Jeans instability anal- 
ysis of Sec. |V] by taking into account the cosmological 
constant A = 310^^^cm^^. This can be done in the 
Newtonian framework by replacing the Poisson equation 
by 



A$ = AnGp 



Ac2 



(CI) 



If we consider an infinite homogeneous distribution of 
matter that is solution of Eq. (|C1[) . we find that the 
density is given by 

A „2 

(C2) 



8ttG' 



Let us assume that these particles are bosons forming 
Bose-Einstein condensates. We restrict ourselves to the 
non-interacting case (a = 0) although more general sit- 
uations could be contemplated as well. The Jeans anal- 
ysis of Sec. |Vl combined with the expression (jC2p of 
the density, shows that the gravitational collapse of these 
self-gravitating BECs is expected to form structures with 
typical radius and mass [see Eqs. ()128|) and ()129|) ] given 
by 



Km?c? 



1/4 



M ' 



(C3) 



Now, using the expression of the cosmological constant 
in terms of classic al fu ndamental constants proposed by 
Bohmer & Harko fl3^: 



A 



^12 ■ 



(C4) 



and introducing the Planck mass, the Planck length and 
the fine structure constant a — /hc= 1/137, we obtain 



-^/p, — (^) Mp. 



12 3/2 



Q,3/2 ^ 



(C5) 



It is interesting to note that these length and mass scales 
can be expressed essentially in terms of fundamental con- 
stants. 



Appendix D: An approximate equation of state for 
relativistic BECs 



Colpi et al. [66[ model a boson star by a scalar field 
with a ^A|(?f)|'' interaction described by the Klein-Gordon- 
Einstein equations. In the Thomas-Fermi limit, they find 
that the scalar field becomes equivalent to a fluid with 
an equation of state 



UK 



1/2 



(Dl) 



where K = Xh'^/Am'^c. For p ^ (small densities), 
we recover the polytropic equation of state p = Kp^ of 
a non-relativistic BEC with short-range interactions de- 
scribed by the Gross-Pitaevskii equation (see Eqs. (fT6)) 
and (|B9p ). For p — )■ -|-oo (high densities), we obtain an 
equation of state p = pc^ I'i like in the core of neutron 
stars. These asymptotic limits were not explicitly given 
in [6^. The study of BECs described by the Tolman- 
Oppenheimer-Volkoff equations with the equation of 
state (jDl|) is considered in |l3l| . It is found that there ex- 
ists a maximum mass M^ax = 0.06136 \f\M^lrr? (very 
close to the result obtained in [g^ by solving the Klein- 
Gordon-Einstein equations) corresponding to a minimum 
radius Rmin ~ 0.3836 ■\/A(A/p/77i) Ac and a maximum 
central density (po)ma2; = 1.592 to'*c'^/A?i'^. 



Appendix E: Conservation of the energy 

The energy associated with the GPP system ([5])- ([7]), or 
equivalently with the quantum barotropic Euler-Poisson 
system ([TO]), ([14]) and © is given by Eq. (HH). Let us 
explicitly show that it is conserved. Using Eq. (f58| . we 
have 



5ec = y —bpdx^ J pu-Sudr. (El) 
Using Eqs. <^ and we find that 

r 

= I ^5pdv = - f QSpdv. (E2) 

2m^ J ^ m J 

According to Eq. (|6l]) . we have 



6U = J H'{p)Spdr^ j h{p)Spdr. (E3) 
Finally, for a symmetric binary potential of interaction 
SW^ [ ^Spdr. (E4) 
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Taking the time derivative of Etot, using the previous 
relations, inserting the hydrodynamic equations (jlOp and 
(fT3| . and integrating by parts, we easily obtain Etot = 0. 
Note that this result remains valid if u is not a potential 
flow since u • (u x (V x u)) = 0. 



Appendix F: Virial theorem 

In this Appendix, we establish the time-dependent 
virial theorem associated with the quantum barotropic 
Euler-Poisson system dTU]), ([HI) and Q. For the sake of 
generality, we derive it in d dimensions. Taking the time 
derivative of the moment of inertia (1721) and using the 
continuity equation (jlOp . we obtain after an integration 
by parts 



1 = 2 / pr ■ udr. 



(Fl) 



With the aid of the continuity equation ([TU|) . the quan- 
tum barotropic Euler equation (fT4)) can be rewritten 



d_ 
dt 



(pu) + V(pu (g) u) 



-Vp - /9V$ - — VQ. 



(F2) 



Taking the time derivative of Eq. (jFip , substituting Eq. 
()F2p . and using integrations by parts, we obtain the time- 
dependent virial theorem 



1 



/ = 2(6^ + eQ)+d pdr + W, 



(F3) 



For a steady state, / = and u = 0, we obtain the 
equilibrium virial theorem 



2eQ + d J pdr + Wu = 0. 



(F4) 



To obtain these expressions, we have used the following 
identities. First, we have introduced the virial of the 



gravitational force 



Wu 



pr ■ V$ dr. 



(F5) 



It can be shown that Wu = (rf — 2)W \i d ^ 2 and 
Wu = —GM'^/2 if d = 2, where W is the gravitational 
(potential) energy |l33l |. On the other hand, it can be 
established after a few manipulations (using essentially 
integrations by part) that 



/ —r-\/Qdr = 2QQ. 
J m 



(F6) 



Indeed, using two successive integrations by parts, we 
have 

d j {V^fdr + 2 j V{r-V^)-V^dr. (F7) 

Using the convention of summation over repeated indices, 
the last integral can be rewritten 

j di{xjdj^)d,y^dr 
= J Xj^^y/p^i^Jy/pdr + j {V dr 
= \j x,d,{V^pfdr + j [V^pfdr 

= j{V^pfdr. (F8) 

Substituting the last term of this expression in Eq. (|F7P , 
we obtain 

y" pr • V {-^^ " ^ / ^^^^^ ^^^^ 

Combining Eqs. (|F9| and (|60l), we obtain Eq. ([F6|) . 



[1] E. Fermi, Rend. Acc. Lincei 3, 145 (1926) 
[2] P. Dirac, Proc. R. Soc. A 112, 661 (1926) 
[3] S.N. Bose, Z. Phys. 26, 178 (1924) 

[4] A. Einstein, Sitzber. Kgl. Preuss. Akad. Wiss. 1, 3 
(1925) 

[5] S. Chandrasekhar, Stellar structure (University of 

Chicago Press, Chicago, 1939) 
[6] S.L. Shapiro, S.A. Teukolsky, Black lioles, white dwarfs 

and neutron stars (Wiley, New York, 1983) 
[7] A.R. Liddle, M.S. Madsen, Int. J. Mod. Phys. D 1, 101 

(1992) 

[8] T.D. Lee, Y. Pang, Phys. Rep. 221, 251 (1992) 
[9] P. Jetzer, Phys. Rep. 220, 163 (1992) 
[10] F.E. Schunck, E.W. Mielke, Class. Quantum Grav. 20, 



R301 (2003) 

[11] R.H. Fowler, Mon. Not. R. Astron. Soc. 87, 114 (1926) 
[12] W.S. Adams, Pub. Astron. Soc. Pac. 27, 236 (1915) 
[13] A.S. Eddington, The Internal Constitution of the Stars 

(Cambridge University, Cambridge, England, 1926) 
[14] W. PauU, Z. Phys. 31, 765 (1925) 
[15] E.C. Stoner, Phil. Mag. 48, 719 (1924) 
[16] E.A. Milne, Mon. Not. R. Astron. Soc. 91, 4 (1930) 
[17] S. Chandrasekhar, Philos. Mag. 70, 592 (1931a) 
[18] L.H. Thomas, Proc. Camb. Phil. Soc. 23, 542 (1927) 
[19] E. Fermi, Rend. Accad. Naz. Lincei 6, 602 (1927); Zeit. 

Phys. 48, 73 (1928) 
[20] S. Chandrasekhar, Astrophys. J. 74, 81 (1931b) 
[21] S. Chandrasekhar, Mon. Not. R. Astron. Soc. 95, 207 



25 



(1935) 

[22] W. Anderson, Z. Phys. 56, 851 (1929) 
[23] E.G. Stoner, Philos. Mag. 9, 944 (1930) 
[24] L. Landau, Phys. Z. Sowjetunion 1, 285 (1932) 
[25] E. Blackman, Nature 440, 148 (2006) 
[26] P.H. Chavanis, Phys. Rev. D 76, 023004 (2007) 
[27] M. Nauenberg, J. Hist. Astron. 39, 297 (2008) 
[28] S.A. Kaplan, Naukovy Zapiski 15, 109 (1949) 
[29] S. Chandrasekhar, R.F. Tooper, Astrophys. J. 139, 
1396 (1964) 

[30] W. Baade, F. Zwicky, Phys. Rev. 45, 138 (1934) 
[31] L.D. Landau, Nature 141, 333 (1938) 
[32] J.R. Oppenheimer, G.M. Volkoff, Phys. Rev. 55, 374 
(1939) 

[33] C.W. Misner, H.S. Zapolsky, Phys. Rev. Lett. 12, 635 
(1964) 

[34] B.K. Harrison, K.S. Thorne B, M. Wakano, J. A. 
Wheeler, Gravitation Theory and Gravitational Gol- 
lapse (Ghicago: University of Ghicago Press, 1965) 
[35] P.H. Ghavanis, Astron. Astrophys. 381, 709 (2002) 
[36] S. Ghandrasekhar, Rev. Mod. Phys. 56, 137 (1984) 
[37] P. Hertel, W. Thirring, Gommun. math. Phys. 24, 22 
(1971) 

[38] J. Messer, J. Math. Phys. 22, 2910 (1981); Phys Lett A 

83, 304 (1981) 
[39] P.H. Ghavanis, Phys. Rev. E 65, 056123 (2002) 
[40] P.H. Ghavanis, Int. J. Mod. Phys. B 20, 3113 (2006) 
[41] V.A. Antonov, Vest. Leningr. Gos. Univ. 7, 135 (1962); 
D. Lynden-Bell, R. Wood, Mon. Not. R. Astron. Soc. 
138, 495 (1968) 
[42] M. Kiesshng, J. Stat. Phys. 55, 203 (1989); P.H. Gha- 
vanis, Astron. Astrophys. 381, 340 (2002) 
[43] M.A. Markov, Phys. Lett. 10, 122 (1964) 
[44] R. Fabbri, R. Jantzen, R. RufSni, Astron. Astrophys. 

114, 219 (1982) 
[45] R. Ruffini, L. Stella, Astron. Astrophys. 119, 35 (1983) 
[46] J.G. Gao, M. Merafina, R. Ruffini, Astron. Astrophys. 
235, 1 (1990) 

[47] N. Bilic, R.D. Viollier, Phys. Lett. B 408, 75 (1997) 
[48] N. Bilic, R.D. Viollier, Eur. Phys. J. G 11, 173 (1999) 
[49] N. Bilic, R.J. Lindebaum, G.B. Tupper, R.D. Viollier, 

Phys. Lett. B 515, 105 (2001) 
[50] N. Bilic, G.B. Tupper, R.D. Viollier, Lect. Notes Phys. 

616, 24 (2003) 

[51] D. Lynden-BeU, Mon. Not. R. Astron. Soc. 136, 101 
(1967) 

[52] A. KuU, R.A. Treumann, H. Bohringer, Astrophys. J. 
466, LI (1996) 

[53] P.H. Ghavanis, J. Sommeria, Mon. Not. R. Astron. Soc. 

296, 569 (1998) 
[54] D.J. Kaup, Phys. Rev. 172, 1331 (1968) 
[55] R. Ruffini, S. Bonazzola, Phys. Rev. 187, 1767 (1969) 
[56] J. A. Wheeler, Phys. Rev. 97, 511 (1955) 
[57] G. Ingrosso, M. Merafina, R. Ruffini, Nuovo Gimento 

105, 977 (1990) 
[58] N. Bific, H. Nikolic, Nucl. Phys. B 590, 575 (2000) 
[59] W. Thirring, Phys. Lett. B 127, 27 (1983) 
[60] J.D. Breit, S. Gupta, A. Zaks, Phys. Lett. B 140, 329 

(1984) 

[61] E. Seidel, W.M. Suen, Phys. Rev. D 42, 384 (1990) 
[62] T.D. Lee, Y. Pang, Nucl. Phys. B 315, 477 (1989) 
[63] M. Gleiser, R. Watkins, Nucl. Phys. B 319, 733 (1989) 
[64] F.V. Kusmartsev, E.W. Mielke, F.E. Schunck, Phys. 
Lett. A 157, 465 (1991) 



[65] E.W. Mielke, F.E. Schunck, Nuc. Phys. B 564, 185 
(2000) 

[66] M. Golpi, S.L. Shapiro, L Wasserman, Phys. Rev. Lett. 
57, 2485 (1986) 

[67] F.E. Schunck, A.R. Liddle, Black Holes: Theory and 
Observation, Proceedings of the 179th W. E. Heraeus 
Seminar. Edited by Friedrich W. Hehl, Glaus Kiefer, 
and Ralph J. K. Metzler. Springer (1998), P. 285 

[68] J.R. Primack, D. Seckel, B. Sadoulet, Annu. Rev. Nucl. 
Part. Sci. 38, 751 (1988) 

[69] J.M. Overduin, P.S. Wesson, Phys. Rep. 402, 267 (2004) 

[70] M.S. Madsen, A.R. Liddle, Phys. Lett. B 251, 507 
(1990) 

[71] M.Yu. Khlopov, B.A. Malomed, Ya.B. Zeldovich, Mon. 

Not. R. astr. Soc. 215, 575 (1985) 
[72] M. Bianchi, D. Grasso, R. Ruffini, Astron. Astrophys. 

231, 301 (1990) 
[73] E. Seidel, W.M. Suen, Phys. Rev. Lett. 72, 2516 (1994) 
[74] G.D. van Albada, Mon. Not. R. Astron. Soc. 201, 939 

(1982) 

[75] P.J.E. Peebles, B. Ratra, Rev. Mod. Phys. 75, 559 
(2003) 

[76] J.F. Navarro, G.S. Frenk, S.D.M. White, Mon. Not. R. 

astr. Soc. 462, 563 (1996) 
[77] A. Burkert, Astrophys. J. 447, L25 (1995); G. Firmani 

et al, Mon. Not. R. astr. Soc. 321, 713 (2001); W.J.G. 

De Blok et al, Mon. Not. R. astr. Soc. 340, 657 (2003); 

G. Gentile et al, Astrophys. J. 634, L145 (2005); G. 

Battaglia et al, Astrophys. J. 681, L13 (2008); S.H. 

Oh et al, arXiv:1011.0899 
[78] G. Kauffmann, S.D.M. White, B. Guiderdoni, Mon. 

Not. R. astr. Soc. 264, 201 (1993); A. Klypin, A.V. 

Kravtsov, O. Valenzuela, Astrophys. J. 522, 82 (1999); 

M. Kamionkowski, A.R. Liddle, Phys. Rev. Lett. 84, 

4525 (2000) 

[79] M.R. Baldeschi, G.B. Gelmini, R. Ruflini, Phys. Lett. 

B 122, 221 (1983) 
[80] S.J. Sin, Phys. Rev. D 50, 3650 (1994) 
[81] R. FerreU, M. Gleiser, Phys. Rev. D 40, 2524 (1989) 
[82] J. Balakrishna, E. Seidel, W.M. Suen, Phys. Rev. D 58, 

104004 (1998) 

[83] F.E. Schunck, astro-ph/9802258' 

[84] T. Matos, F.S. Guzman, F. Astron. Nachr. 320, 97 
(1999) 

[85] F.S. Guzman, T. Matos, Glass. Quantum Grav. 17, L9 
(2000) 

[86] T. Matos, L.A. Ureiia-Lopez, Phys. Rev. D 63, 063506 
(2001) 

[87] M. Alcubierre, F.S. Guzman, T. Matos, D. NMez, L.A. 

Ureiia-Lopez, P. Wiederhold, Glass. Quantum. Grav. 

19, 5017 (2002) 
[88] A. Arbey, J. Lesgourgues, P. Salati, Phys. Rev. D 64, 

123528 (2001) 

[89] M.P. Silverman, R.L. Mallett, Gen. Rel. Grav. 34, 633 
(2002) 

[90] W. Hu, R. Barkana, A. Gruzinov, Phys. Rev. Lett. 85, 
1158 (2000) 

[91] J. Lee, L Koh, Phys. Rev. D 53, 2236 (1996) 
[92] P.J.E. Peebles, Astrophys. J. 534, L127 (2000) 
[93] J. Goodman, New Astronomy 5, 103 (2000) 
[94] A. Arbey, J. Lesgourgues, P. Salati, Phys. Rev. D 68, 
023511 (2003) 

[95] G.G. Bohmer, T. Harko, J. Gosmol. Astropart. Phys. 06 
(2007) 025 



26 



[96; 
[97; 
[98; 
[99; 
[loo; 

101 
102 
103' 
104 
105 
106^ 
107^ 

108; 

109 
110' 

111 
112; 

113 
114 

115 

lie; 

117 

118; 

119' 
120' 

121 
122^ 

123^ 
124 
125 
126; 

127 
128^ 



129' 
130^ 
131 
132 
133^ 
134 



F. Dalfovo, S. Giorgini, L.P. Pitaevskii, S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999) 

A. Bernal, T. Matos, D. Niinez, Rev. Mex. Astron. As- 
trofis. 44, 149 (2008) 

T. Fukuyama, M. Morikawa, T. Tatekawa, JCAP 06, 
033 (2008) 

T. Matos, A. Vazquez-Gonzalez, J. Magana, Mon. Not. 
R. Astron. Soc. 393, 1359 (2009) 

T. Fukuyama, M. Morikawa, Phys. Rev. D 80, 063520 
(2009) 

J.W. Lee, Phys. Lett. B 681, 118 (2009) 
J.W. Lee, S. Lim, JCAP 01, 007 (2010) 

F. Briscese, Phys. Lett. B 696, 315 (2011) 
A. Suarez, T. Matos, ,arXiv:1101.4039] 

P.H. Ghavanis, L. Delfini [arXiv:1103.20"54] (Paper II) 
P.H. Ghavanis, |arXiv: 1103.2698] 

E.P. Gross, Ann. of Phys. 4, 57 (1958); Nuovo Gimento 
20, 454 (1961); J. Math. Phys. 4, 195 (1963) 
L.P. Pitaevskii, Sov. Phys. JETP 9, 830 (1959); ibid 13, 
451 (1961) 

E. Madelung, Zeit. F. Phys. 40, 322 (1927) 

J. Binney, S. Tremaine, Galactic Dynamics (Princeton, 

NJ: Princeton University Press, 1987) 

D. Bohm, Phys. Rev. 85, 166 (1952) 

M. Membrado, A.F. Pacheco, J. Sanudo, Phys. Rev. A 
39, 4207 (1989) 

E. H. Lieb, H.T. Yau, Astrophys. J. 323, 140 (1987) 

E. H. Lieb, R. Seiringer, J. Yngvason, Phys. Rev. A 61, 
043602 (2000) 

G. Manfredi, F. Haas, Phys. Rev. B 64, 075316 (2001) 

F. Haas, G. Manfredi, M. Feix, Phys. Rev. E 62, 2763 
(2000) 

P.H. Ghavanis, Phys. Rev. E 69, 066126 (2004) 

E.B. Kolomeisky, T.J. Newman, J. P. Straley, X. Qi, 

Phys. Rev. Lett. 85, 1146 (2000) 

M. Girardeau, J. Math. Phys. 1, 516 (1960) 

D. D. Holm, J.E. Marsden, T. Ratiu, A. Weinstein, 
Phys. Rep. 123, 1 (1985) 

J. Katz, Found. Phys. 33, 223 (2003) 

A.S. Eddington, Mon. Not. Astron. Soc. 79, 2 (1918) 

P. Ledoux, G.L. Pekeris, Astrophys. J. 94, 124 (1941) 

E. G. Stoner, Philos. Mag. 7, 63 (1929) 

M. Nauenberg, Astrophys. J. 175, 417 (1972) 

G. Sulem, P.L. Sulem, The Nonlinear Schrodinger Equa- 
tion (Springer, 1999) 

P. Sikivie, Q. Yang, Phys. Rev. Lett. 103, 111301 (2009) 
M. Kie sshng, Adv. Appl. Math. 31, 132 (2003); see also 
e-print |arXiv:astro-ph/9910247 
N.N. BogoliubOTTJ. Phys. U.S.S.R. 11, 23 (1947) 
M. Milgrom, Astrophys. J . 270, 365 (1983) 
P.H. 



Ghavanis, T. Harko, |arXiv:ir08. 3986 
G.G. Bohmer, T. Harko, Found. Phys. 38, 216 (2008) 
P.H. Ghavanis, G. Sire, Phys. Rev. E 73, 066103 (2006) 
The observations of Adams (1915) [T^] had led to the 
discovery that white dwarf stars such as Sirius B were 
very compact objects. As noted by Eddington (1926) 
[l3l |: "We have a star of mass about equal to the sun, 
and of radius much less than Uranus." These results 
could not be explained by classical physics. 
[135] As noted by Fermi (1926) H, the Pauli (1925}_[l3] prin- 
ciple was also formulated by Stoner (1924) [l5l |. 
[136] This analogy takes even more sense when one realizes 
that the polytropic model of white dwarf stars that re- 
sults from these considerations UM, uM is mathemati- 



[137] 



[138] 



cally similar to the Thomas-Fermi model of heavy atoms 
developed at the same period [3, [l^ ■ 
For a brief history of the discovery of the limiting mass 
of relativistic white dwarf stars, and the important con- 
tributions of Anderson (1929) Stoner (1930) [U 
and Landau (1932) [S^], some of them having largely 
been forgotten, see Refs. [25l - [27l |. 

As we approach the maximum mass, general relativ- 
ity effects must be taken into account. They induce a 
dynamical instability when the radius becomes smaller 



[139] 



[140] 
[141] 
[142] 



[143] 



[144] 



[145] 



than 1.0267 lO'^'km 



[29[ which is the order of mag- 



nitude of the white dwarfs {R ~ 510 km). The corre- 



1.4176M( 







sponding mass and mean density are M 
and p = 6. 220 lO^gcm"^. 
Two years after Ghadwick's discovery of the neutron in 
1932, Baade & Zwicky (1934) ^ proposed the idea of 
neutron stars and suggested that they would be formed 
in supernova explosions. The possibility of cold, dense 
stars composed principally of neutrons was also contem- 
plated by Landau (1938) [H. 

Self-gravitating degenerate neutrino stars were initially 
suggested as a model of quasars [4^ . 
The case of self-gravitating bosons at nonzero tempera- 
ture has been studied in [st! . [s^ . 

The condensation of bosons of mass m takes place 
provided that the de Broglie wavelength X^b = 
\/2TTh^ /mkBT exceeds the mean separation n^^^'^. 
The critical condensation temperature is fcsTc = 
2nh^n^^^ /{(^{3/2))'^^^m. If we assume an adiabatic cos- 
mological expansion of the universe, the temperature 
has the same dependance T oc n^^^ on the number 
density of the particles. This implies that Bose-Einstein 
condensation occurs if the mass of the particles satisfies 
the condition m < 1.87 eV/c^ [9i.[l33|. 
We note that for a BEG at T = the pressure arising in 
the hydrodynamic equation H14|) has a meaning different 
from the pressure of a normal fluid at flnite temperature. 
It arises directly from the short-range interactions be- 
tween particles encapsulated in the effective potential 
h{p) and is not due to thermal motion (since T = 0). 
As a result, the pressure p can be negative\ This is the 
case, in particular, for a BEG described by the equation 
of state p6|) when the scattering length a is negative 

For boson stars without self-interaction, the Hartree 
(mean-field) approximation, leading to the Schrodinger- 
Poisson system, is exact for A'' — >■ -l-oo (see [ll3t in 
the relativistic case). For self-coupled boson stars with 
fixed a > 0, according to Eq. (|50|) . the TF approxima- 
tion becomes exact for N — -foo. If we want to take 
into account the quantum potential, we have to con- 
sider a nontrivial limit in which A'^ — >■ -foo with aN'^ 
fixed (this differs from the scali ng N — >■ +00 with aN 
fixed considered by Lieb et al. [ll4l ] in the case where 
the gravitational interaction is replaced by an external 
trapping potential). The Hartree (mean-field) approxi- 
mation, leading to the Gross-Pitaevskii-Poisson system, 
becomes exact in that limit. 

Such a description assumes that the fermions have the 
same probability distributions [115| ]. More fundamen- 
tally, the fermions should be described by a mixture 
of A'^ pure states, each of them having a wave function 
^pi obeying the mean-field Schrodinger equation with- 
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out n onlinearity (except the one due to the interaction) 

[nil. 

[146] In the absence of gravitational forces, BECs with a < 

are always unstable. 
[147] The true mass and size of the structures depends on the 

evolution of the system in the nonlinear regime. 
[148] We have seen, however, that spatially inhomogeneous 

distributions with a < are always unstable in the ab- 



sence of gravitational forces. Therefore, the nonlinear 
development of the instability does not lead to equilib- 
rium distributions in that case. We must therefore be 
careful in interpreting the formation of structures di- 
rectly from the linear instability of a spatially homoge- 
neous system. 



